48 | 
VOLUME 81 PAPER No. 748 


ROCEEDINGS 


| AMERICAN SOCIETY 

OF 
| | CIVIL ENGINEERS 
Tt JULY, 1955 


| 
Bll DISCUSSION OF 
Tl PROCEEDINGS PAPERS 
| 1 365, 382, 467, 565, 668 
HYDRAULICS DIVISION 
Copyright 1955 by the American Society oF Civit ENGINEERS 
Printed in the United States of America 
Headquarters of the Society 
33 W. 39th St. 
New York 18, N. Y. 
Tl PRICE $0.50 PER COPY 


SOCIETY OF 
gt 1982 
| 
| 


Current discussion of papers sponsored by the Hydraulics Divi- 
sion is presented as follows: 


Number Page 


365 Effect of Well Screens on Flow into Wells, by 
Jack S. Petersen, Carl Rohwer, and M. L. Albertson. 
(December, 1953. Prior discussion: None. 
Discussion closed) 


Nazir, Mohammad and Ahmad, Nazir 
Rorabaugh, Matthewl. ........ 

Petersen, Jack S., Rohwer, Carl, and Albertson, 


Hydraulic Model Studies of Martin Dam Draft 

Tubes, by Carl E. Kindsvater and R. R. Randolph, Jr. 
(January, 1954. Prior discussion: None. 

Discussion closed) 


Kindsvater, Carl E., and Randolph, R. R., Jr. 


Tranquil Flow Through Open-Channel Constrictions, 
by Carl E. Kindsvater and Rolland W. Carter. 


(August, 1954. Prior discussion: None. Discussion 
closed) 


Laursen, Emmett M., and Toch, Arthur... 
Kindsvater, Carl E., and Carter, Rolland W. 


The Present Status of Research on Sediment 
Transport, by Ning Chien. (December, 1954. 
Prior discussion: 719. Discussion closed) 


Mechanics of Streams with Movable Beds of Fine 
Sand, by Norman H. Brooks. (April, 1955. Prior 
discussion: None. Discussion open until 

August 1, 1955) 


Reprints from this publication may be made on condition that the full 
title of paper, name of author, page reference (or paper number), 
and date of publication by the Society are given. 


The Society is not responsible for any statement made or opinion 
expressed in its publications. 


This paper was published at 1745 S. State Street, Ann Arbor, Mich., 
by the American Society of Civil Engineers. Editorial and General 
Offices are at 33 West Thirty-ninth Street, New York 18, N.Y. 


“eee 5 
8 
10 
“eee 12 
ie ooe 
2 
4 
27 
31 
coos 
565 
... 39 
668 
7 


Discussion of 
"EFFECT OF WELL SCREENS ON FLOW INTO WELLS" 


by Jack S, Petersen, Carl Rohwer, and M, L, Albertson 
(Proc, Paper 365) 


Wen-Hsrune Li," A. M. ASCE.—The authors have conducted a system- 
atic investigation. A simple dimensionless quantity (C L)/D was found for 
determining the hydraulic properties of well screens. 

The equation of continuity, Eq. 5a, is incorrect as the discharge Q in the 
well changes along the direction of flow. The energy equation, Eq. 5b, is also 
incorrect because, even when frictional loss can be neglected, energy is lost 
as a result of the impact of the incoming water on the water flowing in the well. 
That there is a loss of energy caused by impact can be seen from Eq. 12. 
According to the definition of h,, and Eq. 12, 
Q? 


— > > 29 A? + hy. = (Ah,. — Ah’,.) + (ho — Ahy,). . (18) 


in which h, is the piezometric head just outside the screen. The total head 
varies with Ah,, along the flow in the well. Fortunately, Eqs. 5a and 5b 
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Fig. 19.—SckEEN IN AN AQUIFER 


Subscript d Refers to This Level 


were not used in the investigation. With the assumption that the piezometrie 
head h, was constant, Eq. 15 was derived; Eq. 15 agrees with observations on 
screens surrounded by free water (Figs. 2 and 11). The validity of the mo- 
mentum equation, Eq. 5c, and the orifice equation, Eq. 6, is therefore well 
demonstrated. However, one might question the validity of the result for 
screens surrounded by soil in an aquifer. By assuming a constant value 
for h,, the authors have obtained a solution which shows that the discharge 
through the screen is concentrated near the discharge end of the well. If the 
screen were put in an aquifer, such concentration of flow would necessitate a 
decrease of piezometric head h, towurd the discharging end, as shown in Fig. 19. 
The magnitude of the variation of h, depends on the permeability of the soil 
among other factors. The higher the permeability of the soil, the less the 


# Asst. Prof. of Civ. Eng., The Johns Hopkins Univ., Baltimore, Md. 
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variation of h, along the screen. In the tests reported by the authors, the 
free water surrounding the screen is equivalent to a soil with infinite perme- 
ability; therefore, h, is sensibly constant as assumed. When the screen is 
surrounded by soil of lower permeability, h, cannot be assumed to be constant. 

In the analysis of flow into a well, one must account for the permeability 
of the soil in the aquifer. If A is the piezometric head at any point in the 
aquifer, and C, is the coefficient of permeability of the soil (as defined by the 
equation for seepage in isotropic, porous materials), 


Oh 
dh 
q: -C, Oz 
in which q, and q, are the volumes of flow per unit time per unit cross section 
of soil in the r-direction and z-direction, respectively (Fig. 19). Both C, and 
q have the dimensions of length per unit time. As a result of the continuity 
of flow in the aquifer, 
Oh Oh | Hh _ 
or? 


must be satisfied for axial-symmetric flow." To solve Eq. 20, the conditions 
at the boundaries of the aquifer must be known. For the case of a screen 
fully penetrating a confined aquifer (Fig. 19), dh/dz must be zero at the top 
and bottom of the aquifer as g, must vanish there. At great distances from 
the well, dh/dz must also be zero as the flow there is radial. 

The condition at the remaining boundary (just outside the screen) depends 
on the flow in the well. If ha and Q, are h,, and Q, respectively, at the dis- 
charging end of the screen (z = L), the momentum equation, according to 
Eq. 8a, results in 


in which the subscript d refers to the discharge end of the screen. From a 
consideration of the continuity of flow across the screen, 


C, (2) 
Or]. 


in which D, is the outside diameter of the screen, (0h/dr), is 0h/dr atr = D,/2 
and the subscript o refers to the outside surface of the screen. From a con- 
sideration of the continuity of flow in the well, the discharge Q at any point z is 


L L 
dQ, _ oh 


“The Flow of Homogeneous Fluids Through Porous Media,’ by Morris Muskat, McGraw-Hill 
Book C ‘o., Inc., New York, N. Y., 1937, p. 141. 


= 
‘ 
Q? Q? 
A*g 


Substituting the value of Q from Eq. 23 into Eq. 21 results in 


Qs D, C; ah (x D, C,)? / ah 


To relate h,, to the variable h, along the outside surface of the screen, the 
use of the orifice equation, Eq. 6, leads to 


_C. 


Substituting Eqs. 14d, 22, and 24 into Eq. 25 and simplifying, 


| ( 
| (8 D, (a8 
x| ff (3 —~ ), = 


Eq. 26 must be satisfied by the values of h and dh/dr from Eq. 20 just outside 
the screen. The solution of Eq. 20 satisfying the boundary condition specified 
by Eq. 26 can be obtained by relaxation methods."* However, some char- 
acteristics of the flow can be det2rmined without actual numerical computation. 

If hog is the piezometric head in the soil at the discharge end just outside 
the screen, and 6 = hoa — ha (that is, Ah,, at the discharging end), one must 
let@ = h/5,R = 1r/L,and Z = z/L in order to study the problem in dimension- 
less form. Eq. 20 can therefore be reduced to 


1 
dR? " ROR 


Eq. 27 must be solved to satisfy the boundary condition at the screen as 
specified by Eq. 26. Eq. 26 can be rewritten as 


A’g 


1 / a : D (36) 


A’g 


in which the subscript od refers to the discharge end just outside the screen. 
It can be seen that, given the geometrical shape of the aquifer (that is, D,/L) 
and given the coefficients in Eq. 28, a unique solution of Eq. 28 can be obtained. 
That is, the loss coefficient equals 


“ “Relaxation Methods in Theoretical Physics,” by R. V. Southwell, Oxford University Press, Oxford 
England, 1946. 
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The actual value of ¢.4 is immaterial as it represents a constant head for the 
entire aquifer. In order to have some insight into the effect of C, on the loss 
coefficient in Eq. 29, let 


Tee eee eee 


in which m is the ratio of the rate of flow through the screen at any point to 
the mean rate of flow through the screen. With this coefficient m, Eq. 26 


ean be reduced to 
1 2 D 2 
maz +(f maz ) -4(2) m? = 0. .(31) 
CL 


Thus, 


At the discharging end of the screen, Z = 1, h, = hoa, and m = m,. 


A’g 


Thus, the loss coefficient is dependent on oie in which m, is the value of 


m at the discharging end of the screen. The value of m, depends, of course, 
on the independent variables in Eq. 29. In the case of a fully penetrating 
well, the two extreme values of m, are obtained from the case of constant 
h, and the case of uniform discharge through the screen. Both of these cases 
are never realized in the field, but they serve as the limits to the actual cases. 
A sensibly constant h, is obtained with very permeable soils which offer almost 
no resistance to the concentration of the discharge at the discharge end of the 
screen; the authors’ result is an approximate solution for such cases. On the 
other hand, an approximately uniform distribution of discharge through the 
screen is obtained with soils of low permeability which offer high resistance 
to the concentration of the discharge. An approximate solution for such 
cuses can be obtained by setting m, equal to unity. Thus, for small values of 
0,€ 

A*g 


, the loss coefficient can be approximated as 


When Eqs. 15 and 33 are plotted in Fig. 20, it ean be seen that the two curves 
coincide for (C L)/D less than unity. This means that, for all cases in which 
(C L)/D is less than unity, both h, and the discharge through the screen are 
approximately constant. For (C L)/D larger than unity, the curve for a given 
screen may lie anywhere between these two curves. The value of (C L)/D 
for 4 minimum value of the loss coefficient must be larger than 6, and the 
minimum loss coefficient must be less than unity. 
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Additional factors may also influence the value of m., and therefore the 
magnitude of the loss coefficient. For partly penetrating wells, there is a 
tendency for the discharge to concentrate near the bottom of the screen,’ ct 
thus reducing the value of m.. When the screen is surrounded by gravel, the ; 
analysis previously presented must be modified as Eq. 22 is not true when 
there is a redistribution of the discharge through the gravel layer. The actual P 
problem is complicated. The authors’ result gives a loss through the screen 5 8 


ck 
D 


Value of 


A . Value of the Loss Coefficient 


Fic. 20.—Txe Revationsuip Between (CL)/D anv THE Loss 


which is the greatest loss that can possibly occur and is, therefore, of great 
practical value. 

Another point of interest in this subject should also be mentioned. The 
condition of constant A, or uniform distribution of discharge through the 
screen is never realized in nature. Thus, even for fully penetrating wells, the 
flow in the aquifer can never be truly radial. This fact is not generally rec- 
ognized in the study of flow into wells. 


Artuur L. A. M. ASCE.—In the “Synopsis” it is stated: 


“The objective of the investigation was to establish criteria which 
could be used to aid in the selection of well screens to meet the varied 
conditions found throughout the United States.” 


Another interesting statement is contained in the ‘Conclusions’: ‘For the 
criteria to be of greatest practical value to the well driller, a large extension 
of the available coefficients is necessary.”’ It would be of interest to know to 
what extent the laboratory tests should be extended. Also, what value should a 
the well driller place on the test information? q 

In their investigation the authors have anticipated the construction of the 
gravel-packed well. This type of construction necessitates the boring of a 


16“*The Flow of Homogeneous Fluids Through Porous Media,” by Morris Muskat, McGraw-Hill 
Book Co., Inc., New York, N. Y., 1937, p. 269. 


16 Cons. Engr., Berkeley, Calif. : 
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large-diameter well which is protected temporarily by a large casing and then 
by the installation of a smaller casing with perforations. This installation is 
followed by the placing of selected gravel between the casings and finally by 
the removal of the outside casing. 

The reasoning behind the gravel-packed well is correct, but in practice this 
type of construction is unreasonably expensive. In contrast, the well driller 
who is familiar with the substrata in which he is drilling has probably learned 
to construct a well which will satisfy the requirements of the user at a fraction 
of the cost of the made-to-order well with packed gravel. For example, the 
many wells constructed for irrigation purposes are principally located in the 
alluvial fills of vallevs and are not founded in extensive rock formation. There 
must also be present a clay or hardpan stratum on top of the sand and gravel 
through which the water must pass. It is believed that the average well and 
the large-quantity well have clear-water pockets which exist at the top of the 
sand stratum and are bridged by a clay roof. The well, if normal, will have 
part of the perforations entirely unobstructed. As a result of this unobstructed 
condition, there is insufficient outside pressure to cause the water to flow 
through the partly obstructed perforations. 

The well driller uses the factory-perforated casing, or he makes the perfora- 
tions at the proper depth after the casing is in place with a perforating instru- 
ment. By the use of these methods he is able to make a well which is equal 
to, or perhaps superior to, the gravel-packed well at a fraction of the cost. 

In order to explain the water-pocket theory, a description of a well system 
without any perforations will be presented. The well with large perforations 
made after the casing is in place and opposite the assumed pocket is virtually 
a screenless well. 

There are many localities where wells can be constructed without the use 
of screens. Such a location would be one where there are thick clay deposits 
separated by sand and gravel strata of a few feet in thickness. The clay cap 
must be compact enough to bridge the sand if a part of the sand is removed. 
Also, the strata must be below the water level in the well at all times. A well 
that will pump sand to any extent can be assumed to contain a pool in which 
the clear water collects. 

The subject of well development first interested the writer in 1910 when 
he learned that in California ample irrigating water could be obtained from 
bored wells or augered wells. It was also at this time that the deep-well 
turbine became available. This turbine provided the solution to the problem 
of drawing the water from a deep well without providing deep pits for the 
centrifugal pump or for other cumbersome equipment. 

The substrata of the valley floor from Bakersfield to Corning in the San 
Joaquin and Sacramento valleys of California (a distance of approximately 
400 miles) follows this general pattern (clay deposits separated by sand and 
gravel strata). In this area there has been little change in the type of well 
and casing that has been used from 1910 to 1954. Both the slotted section 
made in the shop and the slots made in the field after the well is constructed 
are used. From 1910 until welding equipment became available in the field 
the drillers favored 4-ft-long, double-casing sections. The sections, having 
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slightly different diameters, were telescoped with the joints made at the mid- 
section. The pipe presently used (1954) is the 4-ft section of heavy-gage steel 
with welded collars. 

The strainerless well was patented by the writer on April 28, 1914. The 
idea originated from a clue discovered in a locality where an unusually large 
supply of water could be obtained from a few feet of sand and gravel. In 
these wells the contractor ended the casing some distance above the gravel 
stratum and used no strainer. Elaborating on this idea and applying it to 
other areas where the sands were finer, it was observed that, when the well was 
pumped heavily, sand would be emitted and afterward the water would become 
clear. This action results because the water itself flushes out a percolating 
area until there is a balance between the water velocity and the sand. It was 
also assumed that the percolating space extended radially along the roof of the 
sand stratum from the well to form a sheet of thin depth. 

Therefore, if the casing is extended to the floor of the sand stratum, the 
sand will be moved downward by gravity, and the void area of the well will 
be increased. This increase, however, results from the removal of too much 
sand. If the strainer is opposite the sand, the well is improved by heavily 
perforating the casing near the top of the sand stratum. The strainer adds 
little to the well efficiency. 

After reaching these conclusions concerning the open-bottom well, the writer 
conceived the idea of constructing a well with a 12-in. casing or a 14-in. casing 
which was terminated at the first sand stratum. Then a second string of 10-in. 
casings was extended into the second stratum. The 10-in. casing was then cut 
off near the end of the first stratum—thus leaving an annular space through 
which the water from the first stratum was taken into the casing to join the 
water from the lower stratum. In this way it was possible to extend the open- 
bottom-well principle to several strata in the well. 

It is of interest to note that the writer was a supervising engineer and not 
an experienced well driller when the theory of the open-bottom well was 
deduced. However, not long after the patent was granted, a professional well 
driller sought a license to use the patent. Prior to that date the well driller 
had constructed fifty wells according to this principle, and only when he applied 
for a patent did he learn that one had been previously obtained. 

As a check on the work of this well driller, an agent for a manufacturer 
who furnished most of the well casing for several drillers operating in the area 
observed the wells made by this particular driller. The agent concluded 
that the driller’s wells produced 50% more water (with no sand pumping) than 
did other wells. 

Because of the sight-unseen conditions in a well and the methods of con- 
struction which differ among well drillers, it will always be difficult to obtain 
data—such as the size of the casing, the area of the perforations, and the size 
of the sand grains—which will enable the engineer to predetermine the per- 
formance of a well. 
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Nazir"? anp Nazir Aunmap.'*—A solution to the problem of 
excessive head loss through screens (which are in general use (1954) in the 
Punjab Province of Pakistan) has been presented by the authors. Since the 
inception of the use of ground water, engineers in the Punjab Province have 
used a screen made of }-in.-thick brass sheets in which horizontal slits, 1 in. 
long and of varying widths (from 0.010 in. to 0.014 in.) are cut. These slits 
are sharp on the outside and expand toward the inside at an angle of approxi- 
mately 30°. The distance between adjacent slits is 4 in. In order to maintain 
sufficient strength in the strainer, a }-in. space is left between two sets of 1-in. 
slits. Thus, the open area of these screens is from 6% to 10% of the total 
surface area. 

The prevalent aquifer formation in a part of the Punjab Province is com- 
posed of medium sand; therefore, the coarsest formation has a mean size of 
between 0.30 mm and 0.45 mm. The screens are usually placed in medium 
sand which has a mean diameter of between 0.25 mm and 0.45 mm and an 
effective size of between 0.15 mm and0.25mm. The sand has a high coefficient 
of uniformity—between 1.6 to 2.5—and the usual diameter of the brass screen 


Glass woo! for 
cleaning water 
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pump 
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is between 74 in. and 9 in. A length of approximately 150 ft is considered 
sufficient to pump 2 cu ft per sec. In order to exclude the low-yield region 
and any saline water, a blind pipe with a length of from 40 ft to 100 ft is nor- 
mally installed below the natural surface. 

Several pumping tests have been performed in regions where the water 
table lies from 3 ft to 10 ft below the natural surface. A study has also been 
made of the accompanying lowering of the water table. It has been noted 
that, even at a distance of 2 ft from a well screen, the lowering of the water 
table has seldom exceeded 20% of the pumping level in the well. The lowering 
observed at a distance of 1,000 ft from the well has measured from } ft to 1 ft. 
This small lowering has been attributed to: (1) The formation of free surface; 
(2) the head loss from the formation into the shrouding which envelopes the 
screen; and (3) the head loss between the gravel shrouding and the screen. 


1’ Director, Irrig. Research Inst., Lahore, Pakistan 
Head, Physics Dept., Irrig. Research Inst., Punjab, Lahore, Pakistan. 
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If this lowering is attributed to the concept of a free surface, the existence 
of a long blind pipe and the absence of a discharge face are arguments against 
thisconcept. It can be argued that the small lowering is caused by the greater 
yield of the formation. However, the pumping level in the well contradicts 
this argument. 

Some laboratory experiments have been performed to study the head loss 
between the formation and the gravel and that between the gravel and the 
screen. Assuming the flow in the screen region to be horizontal, the experi- 
mental arrangement shown in Fig. 21 was devised. A supply of air-free 
and filtered water was made available on the inflow side, and on the exit 
side suction up to 20 ft was applied. Small outlets were provided at points 
1, 2, 3, 4, and 5 to record pressure by mercury manometers. The screen 
studied had four different slit sizes—0.014 in., 0.020 in., 0.030 in., and 0.040 
in.—with open areas equal to 7.65%, 11.0%, 13.4%, and 19.3% of the total 
surface area, respectively. The head losses between the gravel and the screen 
and between the sand and the gravel and the observed discharge are shown in 
Table 3. 


TABLE 3.—Tue EFrrects or INCREASE OF SLIT SIZE ON 
DISCHARGE AND Heap Loss 


DiscHarGrE, IN Cvsic 
CENTIMETERS 


Heap Loss, Freer 
PER MINUTE 


Sut Size, 1x Incues 


Suction | —— 


ae. ig 014) 0.020) 0.030) 0.040 0.014 0.020 0.030 0.040 
Open area, in percentage | | 
of total area Between Caused | tween) Caused|Betwee Caused 
sand | sand sand | sand b 
and and | and and y 
| ] grave | screen | gravel | screen gravel screen gravel screen 
| 7.65 | 11.01! 13.4 | 19.3 | | 
5 224 | 238 260 | 266 0.62 0.17 | 058 | 017 | 046 | 009 | 037 | 0.08 
7 302 | 318 | 346 | 354 1.22 | 0.20 0.94 | 0.20 0.76 | O12 | 0.64 0.10 
10 421 43606476 | 485 1.82 | O 23 160 | 0.23 | 140 | 0.16 | 1.26 0.15 
12 | 504 | 519 560) 568 218 | 026 | 200 $025 | 188 | 0.20 1.70 0.18 
15 | 636 | 640 | 692 | 710 2.65 | 0.28 | 245 | 0.28 2.50 0.15 | 2.32 | 0.30 
17 | 720 | 722 | 784 | 808 290 | 0.30 2.70 | 0.30 2.81 | 0.31 2.74 0.30 


A study of Table 3 reveals that the head loss between the screen and the 
gravel is nearly constant, the variation being from 0.1 ft to 0.3 ft with suction 
head varying from 5 ft to 17 ft. The head loss between the formation and the 
gravel is, however, significant, as it is approximately ten times as great as that 
between the gravel and the screen. Similarly, the discharge with a given 
suction head was the highest for the screen with the widest slits. Experiments 
with wider slits and open area greater than 20% are in progress (1954), but 
the preliminary observations have shown that screens with slits wider than 
0.040 in. may not increase the discharge under a given suction head. 

Although there is no question about the mathematical derivation of Eq. 15, 
it does seem that in the authors’ range of application, laminar flow should exist 
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and there should be no movement of particles. Therefore, the effect of viscosity 
cannot be eliminated. Similarly, the statement that the factor (C L)/D should 
always be greater than 6 to reduce the screen loss to a minimum seems 
difficult to comprehend. From the writers’ experience in Pakistan, D does 
not exceed 9.0 in., L can be as great as 150 ft., and the value of C (from Eq. 
14d) can be much greater than 30. Therefore, the factor (C L)/D will 
always be greater than 6, and the loss will be minimum regardless of the water 
pumped out. However, a consideration of this point does not explain the high 
pumping level, high yield, and attendant slight lowering of the water table in 
Pakistan. 


Matruew I. RorasavuGn,’” A. M. ASCE.—The experimental data pre- 
sented by Messrs. Petersen, Rohwer, and Albertson are a valuable contribu- 
tion to the hydraulics of wells. Much speculation and many claims have 
been made concerning the magnitude of head losses through well screens. 
The data presented should do much toward clarifying some of the problems 
involved in screen design and in the selection of well screens. 

The authors are to be complimented on their work considering the difficulties 
inherent in such a study. The economical construction of a suitable model 
which will be hydraulically equivalent to the prototype is particularly difficult. 
Small-scale models of screens are costly items. Other items also add to the 
cost of research and to the sensitivity of the necessary physical observations. 

One of the basic assumptions in the theoretical development is that the 
piezometric héad on the outside of the screen is constant over the entire outside 
surface of the screen. In the case of a gravel envelope, this condition was not 
met, as is implied in Fig. 8. The authors considered this variation to be rela- 
tively small and concluded that the status of the investigation did not warrant 
the use of a correction for the variation. The effects shown in Fig. 8 result, in 
part, from the fact that the screen (2 ft long) does not fully penetrate the satur- 
ated thickness (5 ft). Flow enters the gravel pack over this 5-ft length and 
travels through the gravel to discharge through the 2-ft screen. Flow lines 
converge radially in the horizontal plane and also converge in the vertical plane. 
The head losses in Fig. 8 are taken horizontally through the gravel pack. The 
shape of the curves in Fig. 8 is partly controlled by the position of the screen in 
relation to the boundaries—the bottom of the tank and the water surface in the 
tank. 

The effects of the boundary conditions on flow toward the screen in the 
laboratory structure are related to the data in Figs. 16, 17, and 18. From 
reference to the test apparatus shown in Fig. 6 and to the tests on screens with- 
out a surrounding grave pack, it should be evident that, as A, approaches 
100%, the horizontal component of velocity at the screen varies from approxi- 
mately zero at the water surface to a maximum near the base of the screen. 
This generalization can be visualized easily by considering the flow net that 
would be developed between the overflow rim and the outflow pipe shown in 
Fig. 6. That is, for the condition of no screen or screens having a low friction 
loss, head-distribution curves similar in shape to those in Figs. 16, 17, and 18 
would result. These curves would be based on the flow pattern in the tank 


® Dist. Engr Geological Survey, U s Dept. of the Interior, Ground Water Branch, Louisville, Ky. 
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and would be controlled by the tank boundaries and the position of the inflow 
and outflow pipes. Item 4 in the ‘‘Conclusions,’”’ apparently based on those 
data in Fig. 17, applies to the test apparatus used. However, considering the 
effects of boundary conditions and of the partial penetration introduced by the 
addition of a gravel pack 9 in. thick (as demonstrated in Fig. 8), it is doubtful 
whether this conclusion can be applied to prototype structures. 

The investigation deals with head losses in a very small zone in and near the 
screen. This is only a small part of the problem of total drawdown in a well. 
Before the screen data are used to form conclusions for general well design, it 
would be desirable to consider comprehensively the head distribution in and 
about the pumping well. 

If an artesian gravel bed, 4 ft thick, having a permeability of 2,500 gal per 
day per sq ft and a storage coefficient of 0.003, is fully penetrated by a 4-ft-long, 
3-in.-diameter, A-d-type, well screen, and if pumping is conducted at rate of 
0.125 cu ft per sec, the conditions approximate those shown in Fig. 17. Fora 
pumping time of 1 hr the drawdown in the well, computed by the nonequilib- 
rium formula developed” by C. V. Theis, is approximately 10 ft, if turbulent 
losses in the aquifer are not considered." According to this conclusion (item 
4) the most efficient well screen would have to penetrate only 1 ft of the 4-ft 
aquifer (Fig. 17). Using the Kozeny* equation for partial penetration, the 
drawdown at the well is approximately 15 ft rather than 10 ft. A plot of head 
distribution about a partly penetrating well is given by M. Muskat.” Thus, 
the selection of an efficient well screen for this example (by use of item 4 in the 
‘‘Conclusions’’) results in a crowding of the flow lines outside the well and in- 
efficient use of the aquifer. Although the head loss for the 1-ft screen and the 
4-ft screen will be almost the same when tested in a tank, the drawdown will 
differ by approximately 5 ft when installed in an aquifer of the type used in the 
example. 

The permeability of the aquifer is an important item in examining the vari- 
ous factors influencing well design. In an aquifer having a very high permeabil- 
ity, the aquifer losses will be small, and attention to a modified form of item 4 in 
the ‘“‘Conclusions”’ could be justified by a consideration of the head distribution 
in the aquifer. However, in most of the unconsolidated aquifers tested in the 
United States, the usual operating head loss is measured in feet and in tens of 
feet. For these wells, a screen loss of a few tenths of a foot is small in propor- 
tion to the total loss, and item 4 in the ‘‘Conclusions’”’ would not form an im- 
portant part of the economies in the entire flow system. 

It is noted that the rates of flow used in the studies were, in general, higher 
than those in most producing wells. Rates in the tests, for example, were 
0.125 cu ft per sec through | ft of a 3-in. screen and 1 cu ft per sec through 2 ft 
of a 12-in. screen. With these rates, only a few tenths of a foot of head loss 

2%**The Relation Between the Lowering of the Piezometric Surface and the Rate and Duration of Dis- 
Vale ee oe Ground-Water Storage,"’ by C. V. Theis, Transactions, Am. Geophysical Union, 


"Graphical and Theoretical Analysis of Step-Drawdown Test of Artesian Well,” by M. I. Rora- 
baugh, Proceedings-Separate No. 8362, ASCE, December, 1953. 


2**The Flow of Homogeneous Fluids Through Porous Media,"’ by M. Muskat, J. W. Edwards, Inc., Ann 
Arbor, Mich., 1946, p. 274. 


Ihid., p. 270. 
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was measured. In a field installation, using 20 ft of a 12-in. sereen and pro- 
ducing 500 gal per min, losses through the screen become small when compared 
to other losses. 

This investigation should do much to clear up misunderstandings between 
owners, drilling contractors, and screen manufacturers as to the extent of the 
head losses of screens. The data demonstrate that losses in most commercial 
screens are small, and that the screen manufacturers have done well in designing 
efficient screens. 


Dean F. Peterson, Jr..4 M. ASCE.—A most interesting experimental 
verification of a theoretical hypothesis in fluid mechanics has been presented ; 
experimental results are seldom as satisfactory as in this ease. The authors 
are to be complimented on their application of the basic laws of hydrodynamics 
to the development of this hypothesis and the manner in which they organized 
their experimental program to check it. 

In the tests that omitted a gravel pack, the outside piezometric head was 
essentially constant. This constancy could not be true of a well in a sand 
formation or in a gravel formation. The screen phenomenon would un- 
doubtedly cause the streamlines approaching the well to converge near the 
discharge end of the screen and would thus cause pressures outside the casing 
to be lowest near the discharge end. This would tend to overcome, to some 
extent, the effects of the screen which were reported by the authors. The 
foregoing effect would certainly be secondary, and it is expected that the authors’ 
analysis will still give useful quantitative results. 

Most welcome would be a solution to the problem of the hydraulics of a 
well in an aquifer which would take into account the screen effects suggested 
by the authors. All solutions of the problem of well hydraulics, such as that 
offered by Jules Dupuit,?® are based on the assumption of equal piezometric head 
throughout the length of the well. 

Another interesting problem is the effect of a pump intake located within 
the length of the screen. The velocities inside the screen would converge from 
both directions toward the pump intake. One would thus infer that, by 
proper placement of the pump intake, a “double effective” length of screen 
might be developed—one length above the intake and one below. For the 
most effective operation, it would appear that the pump intake would need 
to be set above the bottom of the screen a distance such that (C L)/D is greater 
than 6 in order that the screen might be unrestricted by pump bowls and 
columns throughout the effective length. In addition, there would be an 
effective distance above the pump intake of (C L)/D equal to 6. For this latter 
‘ase, D would need to be corrected to allow for the area occupied by the pipe 
columns and bowls. Experiments on an effective screen length using actual 
pump installations would be enlightening; perhaps this might be done by the 
use of a small-scale pump and screen. 


* Prof. and Head of Civ. Eng., Colorado Agri. and Mech. College, Fort Collins, Colo. 
* Etudes théoretiques et pratiques sur le mouvement des eaux,"” by Jules Dupuit, Paris, France, 


2d Ed., 1863 
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G. Conen ve Lara.2*—The development of ground-water resources in 
Algeria and Tunisia necessitated a general study of friction losses in well screens 
and in the surrounding aquifer. This study, begun in 1949, resembles in part 
that performed by Mess:s. Petersen, Rohwer, and Albertson. On the whole, 
these results agree with those of the authors, but some of their conclusions ap- 
pear to be incomplete. 

The problem was approached by determining the theoretical laws of con- 
tinuous flow through screens, with and without porous material surrounding 
them. In order to support the theory, air-flow tests (without surrounding 
material) were performed according to the Reynolds similarity. The tests 


ee eee vee 


Fig. 22.—Sewictrcr tar Mover 


with surrounding material were conducted with water on the semicircular model 
shown in Fig. 22. The results of these latter tests have not been published 
(as of 1954). The purpose of these tests is similar to the purpose of those per- 
formed by the authors, although different methods were used. 

Theoretical Computation and Experimental Research for a Screen Surrounded 
Only by a Liquid.—The theoretical computation is similar to that offered?’ by 
P. Arnaud, and the experimental measurements accurately confirm the basic 
assumptions. It is agreed that the experimental determination of a coefficient 
C, (which varies according to the type of commercial screen used) allows one 
to apply the theoretical computation to commercial screens, to which the 
basic assumptions used in the computation do not always apply. 


2¢ Engr., Laboratoire Dauphinois d’' Hydraulique, Grenoble, France 


27 *Pertes de Charge dans les Crépines de Puits de Pompage,” by P. Arnaud, Troisiéme Journées de 
ydraulique, Paris, France, 1954 
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Well Screen Surrounded by Gravel.—During the tests on a circular model, it 
was found (as also by the authors) that the pressure on the exterior wall of the 
screen was not constant. For a louvre-type screen (200 mm in diameter and 
70 cm in length) with high discharges, there resulted curves similar to those 
shown in Fig. 8. It is not known whether the trend of these curves is main- 
tained for small discharges as the resulting pressure differences were too small 
to be measured. 

In the theoretical computation, it is supposed that the pressure on the ex- 
terior wall is constant. The authors have tried to avoid the difficulty by deter- 
mining a coefficient, C,, similar to the one used for commercial screens sur- 
rounded only by liquid. The head loss, Ah, entering in the dimensionless 


number, an |x , was taken as the difference in pressure inated the static 


pressure intake in the interior of the screen and the static pressure intake half 
way along the screen, between the gravel and the screen—that is, where there is 
a minimum static pressure for the given discharge. Based on this definition, the 
authors determined several values for C,, which appears as a function of the 
opening dimensions and size of gravel used for a given type of screen. It is 
implied by the authors that C, is applicable regardless of the voids ratio of the 
surrounding soil and of the thickness of the gravel. This result seems incom- 
plete as a change in the surrounding soil thickness or of its voids ratio will entail 
a variation of Ah—therefore of C,—even for the same discharge using the same 
apparatus. 

Furthermore, a second theoretical computation has been established?’ by 
Mr. Arnaud for a screen surrounded by gravel in which head losses are assumed 
to be laminar. This computation has been verified by tests on the semicircular 
model. The results of this theoretical and experimental work show that the 


2 
minimum length, L, of the screen—beyond which an /F remains constant 


regardless of further increase of the screen length—depends on the permeability, 
K, and the thickness of the surrounding material. The resulting formula 
is 


AH _ {LK vk 


in which AH is the drop in head across the screen; K denotes the permeability 
of the soil; ¢t is the thickness of the surrounding soil; r represents the screen 
radius; h is the total head loss through the screen, gravel, and aquifer; and 
C,, denotes the discharge coefficient for the screen slots. The surrounding 
material also changes the contraction coefficient for the openings and the per- 
centage of open area of the screen. These two changes depend on the relative 
dimensions of screen openings and gravel. The over-all effect of the surround- 
ing gravel on the screen, characterized by the discharge coefficient, C,, (Eq. 
34), was determined by a comparison—at constant discharge—of head losses 
in a screen and filter of given thickness obtained in the two following cases: (1) 
A well screen surrounded by fine wire gauze, with a percentage of open area 
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higher than that of the surrounding material, and (2) an ordinary well screen, 
without wire gauze. The discharge coefficient, C,,, could thus be experiment- 
ally determined as a function of the mean diameter of the surrounding grains. 

The values of C, determined by Messrs. Petersen, Rohwer, and Albertson 
therefore apply only to the conditions on their model—that is, for the filter 
thickness and voids ratio used by them. 


TABLE 4.—ReEvaTionsHip AMONG SCREEN COEFFICIENT, GRAVEL 
SizE, AND SLot WipTH 


Mean diameter Slot width, in Mean value 
Screen coefficient | of gravel, in inches inches Col. 2 + Col. 3 of Col. 4 
(1) (2) (3) (4) (5) 
ty 
Jl 1.2 
0.3 0.23 * 1.23 1.23 
q 03 4 1.2 
0.4 0.36 * 1.92 1.9 
4 0.44 4 1.75 
‘ 
0.5 0.55 2°94 2.75 
0.68 4 2.71 
f 0.24 te 3.84 
0.6 0.55 t 49 47 
0.95 49 


* Values obtained from Fig. 13. 


Value of the Screen Coefficient 


Fro. 23.—VartaTion or Screen Corrrictent with Graver Size ann 


The Variation of the Value of the Screen Coeffiicient.—Fig. 13 represents the 
variation of C, as a function of the size of surrounding gravel for different widths 
of type B screen openings. It appears from Fig. 13 that a geometric similitude 
does not exist and that different values of C, can be found for a given ratio of 
d (diameter of gravel) to s (slot size). Using the authors’ results, it was 
possible to see that all the curves of Fig. 13 can be reduced to a single curve, 
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with the variation of C, considered to be a function of d/s (Table 4 and Fig. 23), 
Geometric similitude is thus obtained for type B screens. It would be inter- 
esting to make the same observations for screens of other types. 

Practical Application.—The practical application of the results obtained in 
this study requires a knowledge of the hvdro-geologic variables entering into 
Eq. 34 and of the permeability of the surrounding aquifer. The composition of 
a natural or artificial filter, intended to prevent sand particles from entering 
the well and to diminish the loss of head in the immediate vicinity of the well, 
complicates the problem. The minimum length that can be given to a screen 
will depend on the thickness of the natural or artificial filter, the hydro-geologic 
conditions, and the permeability of the filter and aquifer. If it were possible 
to determine C, as a function of these variables and to determine the effect of 
C, on C, it would be possible to apply Conclusion 1 of the paper. 

In view of the difficulties that must be overcome, it can be seen how much 
knowledge and intuition are needed to determine accurately the minimum 
length of a well screen. Nevertheless, an examination of the combination of 
screen and soil shows that a great part of the head loss occurs in the immediate 
vicinity of the screen because of the convergence of flow. It is then possible to 
determine an “effective radius of influence,” which can be defined as the radius 
of the zone within which 95% of the total head loss in the combination screen- 
aquifer occurs. Knowing the soil composition and assuming maximum con- 
solidation, it is possible to determine the minimum economical length of the 
screen to be installed in a well. 


Gérarp Tison, Jr.2*.—At the hydraulic laboratory of the University of 


Ghent in Belgium, some investigations have been made concerning the resist- 
ance to a flow of water through a special type of well sereen. The screen was 
made of steel with large elongated holes and was surrounded by an envelope 
(of gravel or sand) 14 mm (0.55 in.) thick. The grains of this sand or gravel 
were bound by plastic matter which was kept to a minimum so that the per- 
meability of the envelope was held as nearly as possible to the natural perme- 
ability. The grain size of the sand used in the test sample ranged from 0.7 
mm to 1.2 mm (from 0.028 in. to 0.047 in.). 

The methods of investigation were slightly different from those used by the 
authors. A well-screen element 0.86 m (2.62 ft) long and with a 0.076-m 
(1.93 in.) inside diameter was placed along the horizontal axis of a steel cylin- 
drical tank the diameter of which was 0.90 m (2.74 ft) and the length of which was 
1.25 m (3.81 ft). This cylindrical tank (Fig. 24) was closed at its two ends by 
steel sheets. A 2-in.-diameter pipe which was fitted into one of the end plates 
of the cylinder was the outlet for the flow through the screen. 

The cylindrical tank was connected with a water supply. The pressure in 
the tank was maintained at a constant magnitude by a valve on the supply 
pipe. Under the effect of that pressure the water flowed through the well 
screen. Outside the well screen there was hydrostatic pressure. Inside the 
well screen the value of the pressure was a maximum at the closed end and a 


minimum at the outlet pipe. 
* Civ. Engr., Hydr. Lab., Univ. of Ghent, Ghent, Belgium. 
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Investigations concerning filtration have shown that movement is laminar 
so long as 
V,d 


in which Vy, the filtration velocity, is 
Q 
( 


R represents the Reynolds number, d is the effective grain size, v denotes the 
kinematic viscosity, and A, is the area of the filtration surface. 


4 


Fic. ExpertMentaL Went SCREEN 
Several of the investigations at the University of Ghent satisfied the condi- ae 


tion represented by Eq. 35. It is known that in this case the discharge is 
proportional to the difference of piezometric head between the inside and the 

outside of the screen. Applying the momentum equation to the flow within 
the pipe between the closed section (0) through which no flow passes and any ‘ 
section (1) of that pipe, | 


A 
As the piezometric head is constant outside the well screen, Eq. 37 can be | 
written as 
Differentiating Eq. 38 with respect to the length L results in 
dQ Ag d(Ah,:) _ Vg __ d(Ahy:) (39) 


However, the increment of discharge dQ into the screen corresponding with an 
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increment of length dL is 


dQ = K™? db 


in which e¢ is the thickness of the gravel envelope. 
Combining Eqs. 39 and 40 results in 


VgeD d(Ah,:) | 


dL = 


Integration of Eq. 41 yields 


in which the constant of integration equals zero. The value of K for the en- 
velope was found to be 0.00525 m per sec. With R equal to 5, it was possible 
to measure Ah,, in certain sections inside the screen well by using a pitot tube. 
The results of these measurements (and a comparison of the measured values 
with those obtained from Eq. 6) are as follows: 


Distance between 
closed end of Value of Ahp: 

drain and section, from Eq. 42, Measured value of 
in centimeters in meters Shpz, in meters 


0.0166 


0.0182 
0.019 

0.0195 
0.0202 
0.0215 
0.0229 
0.0246 
Rs 0.0262 
0.0275 
0.0278 
0.028 


0.018 
0.019 
0.0195 
0.02 
0.0215 
0.023 
0.025 
0.0265 
0.028 
0.03 
0.05 


The previously tabulated discrepancies in the value of Ah,, are negligible 
except for the last sections close to the outlet tube. The discrepancy can be 
explained by the contraction of flow in the vicinity of the outlet. This con- 
traction increases the velocity and decreases the head inside the tube. Assum- 
ing for Ah,. the value 0.0275 m (the last measurement on which the contraction 
has no influence), it is found from Eq. 40 that Q = 0.00148 cu m persec. The 
measured discharge was 0.001495 cu m per sec. 

When R is greater than 5, laboratory tests have shown that the filtration 
movement is no longer laminar. The discharge is then proportional to (A4h,,.)" 
with n having values between 1 and 3. Hence, integration of Eq. 41 becomes 
impossible. 
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Jack S. Petersen,” Cart Ronwer,” anp Maurice L. ALBERTSON,” 
Memsers, ASCE.—Certain interesting ideas worthy of serious consideration 
have been advanced by Mr. Li. Contrary to Mr. Li’s statement, Eqs. 5a and 
5b are correct because they are given simply as the ‘“‘usual forms.”’ Unfor- 
tunately, it was not made clear that, when Eqs. 5a and 5b are used, the sub- 
scripts do not correspond with those in Fig. 1. When applied to the flow 
approaching and passing through the screen opening and to the jet issuing from 
the opening, the subscripts must be changed accordingly. Furthermore, the 
energy which is lost is that of the jet, not that of the longitudinal flow, as im- 
plied by Mr. Li. Therefore, Eqs. 5a and 5b were used in deriving the orifice 
equation, Eq. 6, which was used to obtain Eq. 14a. 

Mr. Li states that ‘“* * * one might question the validity of the result for 
screens surrounded by soil in an aquifer.”” It is true that the higher the visco- 
sity of the fluid and the lower the permeability of the aquifer, the greater the 
variation of the piezometric head along the outside of the screen. However, 
for simplicity in deriving the theoretical equations, this variation was assumed 
to be of secondary importance. The data (shown in Fig. 11) clearly demon- 
strate that this assumption was justified. 

Within the aquifer at a substantial distance from the screen the flow is 
undoubtedly laminar, as it must be if Eqs. 19a and 19 are to be applicable. 
However, it is possible that for certain conditions the Reynolds number within 
the gravel envelope is sufficiently high to cause transition flow or even turbulent 
flow, in which case Eqs. 19a and 19) are not applicable. If the Reynolds num- 
ber is greater than 1.0, the flow is no longer laminar; in computing the Reynolds 
number, D is the grain diameter and V is computed from the bulk area. 

Eq. 29, as derived by Mr. Li, is an excellent contribution to the investiga- 
tion. As in the case of Eq. 16, however, such a derivation is no stronger than 
its weakest assumption. The final test is its comparison with the data. There- 
fore, experimental data are necessary to establish the exact influence of each of 
the parameters in Eq. 29. 

In Eq. 32, m is assumed to be independent of the loss coefficient and the 
screen coefficient, and Eq. 33 further Gepends on m being equal to unity—all of 
which are not necessarily true. The writers have evaluated all the data at their 
disposal (including experiments they have made using a gravel enveiope), and 
in no case is the loss coefficient less than 1.0. Therefore, despite the ingenious 
derivation of Eq. 33 and the desire of the writers to utilize it, it must be con- 
cluded that either the limiting assumptions make it invalid or the range of 
experimental data is not sufficiently wide to include this special case. 

Although the open-bottom wells described by Mr. Collins are used success- 
fully in coping with the problem of sand pumping when suitable clay or rock 
strata cover the aquifer, there are many areas where these strata are not strong 
enough to support the overlying material. Where this condition exists other 
methods of constructing wells have to be adopted. Other types of wells may 


2% Instructor, Univ. of Missouri, Columbia, Mo. 
*® Fort Collins, Colo. 


® Prof. of Civ. Eng. and Head of Fluid Mechanics Research, Colorado Agri. and Mech. College, Fort 
Collins, Colo. 
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sometimes develop large cavities similar to those produced in open-bottom 
wells; this condition is dangerous because collapse of the roof over the cavity 
may ruin the well. One reason for the use of gravel-packed wells is that the 
gravel tends to prevent the formation of these cavities. 

It was not the intention of the writers to recommend the construction of any 
one type of well but, because wells with screens (with or without gravel envel- 
opes) are most common, a study of the hydraulic principles governing the flow of 
water into such wells seemed desirable. 

As to the questions raised by Mr. Collins regarding the tests, it is believed 
that well drillers should take account of the facts established by the laboratory 
experiments. The writers, however, are aware that field confirmation of the 
proposed theory of flow is desirable even though the laboratory tests were in 
remarkably good agreement with the theory. 

The writers are interested in the experiments on the loss of head through 
screens and gravel packs reported by Messrs. Nazir and Ahmad. The report 
on the results of loss-of-head tests on screens with larger percentages of open- 
ings should be most valuable. Although the results reported in Table 3 show, 
in general, that the head losses through the screens decrease as the percentage of 
open area increases, the tests on the screens with larger openings should demon- 
strate that the losses hecome constant. The data on the loss of head through 
the gravel for a constant suction head indicate that the loss decreases as the 
discharge increases for the larger slit openings. It is difficult to understand 
why this decrease should occur unless fine material was being washed out of the 
gravel. 

Although the wells in Pakistan have values of (C L)/D greater than 30, it 
does not necessarily follow that all of the screen in the well is effective in de- 
creasing the loss of head. If the writers’ analysis is valid, only that part of the 
screen which will make (C L)/D = 6 will be utilized (Fig. 17). However, it 
was not their intention to imply that the loss of head would be the same mini- 
mum regardless of the discharge. For each discharge there is a minimum loss 
which an increase in the length of screen or the percentage of openings will not 
reduce, 

Mr. Peterson’s suggestion that a “double effective’ length might be de- 
veloped in a well by properly locating the suction inlet seems reasonable. This 
phase of the problem was not investigated, but the flow from above the suction 
inlet—as well as the flow from below the suction inlet—should follow the cri- 
terion presented by the writers. 

The experiments performed at the University of Ghent, as reported by Mr. 
Tison, take into account the effect of a special gravel envelope on the head 
losses throug’: the screen. Although the equipment and the analysis differed 
from those used by the writers, the head losses follow the pattern reported by 
them, as shown in Figs. 16, 17, and 18. Unfortunately, the percentage ‘of 
openings and the coefficient of discharge are unknown for this type of screen. 
Otherwise it would have been possible to compute the length of screen required 
to reach the limiting value of (C L)/D = 6. 

The writers are of the opinion that, if the (AA,.)-values for distances of 82 
em (32.3 in.) and 86 em (33.8 in.) between the closed end of the drain and the 
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section had been measured, still larger head losses would have been observed. 
It is also believed that the observed losses reported by Mr. Tison may be 
more nearly correct than the computed losses. 

Problems involving many variables can be approached in many ways; 
different investigators may derive formulas of different form for the solution of 
the same problem by using a different approach. Eq. 34 reported by Mr. de 
Lara, provides for the thickness and permeability of the gravel envelope 
whereas Eq. 15 takes account of the gravel envelope in so far as it affects the 
value of C by changing the orifice coefficient and the effective open area of the 
screen. As shown in Figs. 14 and 15, the loss through the screen is independent 
of the gravel envelope except when the percentage of open area in the screen is 
inadequate and the size of the gravel is small. 

An important relationship based on Fig. 13 is shown by Mr. de Lara in 
Table 4 and Fig. 23. Fig. 23 shows that the curves for the different slot 
openings in Fig. 13 can be combined into a single curve by substituting the 
ratio d/s for the gravel size. However, when d/s is less than 1, the gravel 
particles would not be held back effectively by the screen, and an unstable 
condition would result. For this reason the part of the curve for d/s less than 
1 should either be eliminated from Fig. 23 or be shown as a broken line. No 
doubt the screen coefficients for the type C screens could also be plotted 
theoretically on a single curve because of the geometrical similitude of the slots 
in these screens. However, these screens have a relatively large percentage of 
open area, and as a result most of the tests were in the range in which the head 
losses were a minimum—that is, (C L)/D greater than 6. Consequently, 
sufficient data are not available for the analysis because the effective length of 
the screen is unknown. 

The writers recognize that the conclusions stated in the paper require con- 
firmation by field tests. Their primary purpose was to call attention to the 
peculiarities of flow into wells that had not been previously reported. 

As noted by Mr. Rorabaugh, the data are based on tests of screens that do 
not penetrate the aquifer completely. The original plan for the tests included 
the installation of flexible disks in the gravel at the top and bottom of the 
perforated area of the screens so that all the water would approach the screen 
radially and horizontally. This plan was not satisfactory because the disks 
moved as the gravel settled. The disks also produced a plane of different 
permeability at the top and bottom of the screen. For these reasons the disks 
were not used. 

It should be noted that complete penetration of the screen would not have 
eliminated the converging flow lines because the tests showed that most of the 
water entered the screen at the discharging end. This condition tends to 
concentrate the flow lines at the discharge end. The same effect is shown by 
the tests reported by Mr. Tison and by Mr. Li in Fig. 19. 

The head loss through the screen and the gravel pack is a small fraction of 
the total drawdown because most of the head loss occurs in the aquifer. The 
loss in the aquifer depends on the discharge and the characteristics of the 
aquifer and not on the screen or the gravel pack. 
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In Mr. Rorabaugh’s analysis of the drawdown in a well in an artesian aquifer 
under conditions of full and partial penetration it is shown that the drawdown 
decreases as the depth of penetration increases. No doubt this is true under 
certain conditions, but the writers believe that it may not be true if the length 
of screen required to make the value of (C L)/D equal to 6 is less than the depth 
of penetration. So far as is known, these two conditions have never been 
checked by actual tests. Such a test is possible and would certainly yield use- 
ful results. 

The writers appreciate the interest shown by all those who have contributed 
discussions. It is hoped that this work will lead others to investigate the prob- 
lem by conducting field tests wherever conditions are favorable. 


Corrections for Transactions.—Eq. 8b should be changed to — A’gdhy,, 
= d(Q*). Eq. 9a should be changed to Sh,, = a constant —h,,. Eq. 9b 
should be changed to d(Ah,,) = — dh,,. In Eqs. 14a and 146 the radical 
signs in the denominators of the right-hand members should be extended to 
cover the terms, Mh’,,. All parentheses should be removed from Eqs. 15 and 


16. 
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Discussion of 
"HYDRAULIC MODEL STUDIES OF MARTIN DAM DRAFT TUBES" 


by Carl E, Kindsvater and R, R, Randolph, Jr, 
(Proc, Paper 382) 

Ricwarp §. Wooprurr,* A. M. ASCE.—Pyramids were installed in the 

§ Senior Engr., Alabama Power Co., Birmingham, Ala. 
three original draft tubes in September, October, and November, 1952, and 
were of the dimensions recommended on the basis of the authors’ research (Fig. 
11). When the draft tubes of Units 2 and 3 were unwatered, large irregular 
pieces of concrete were found immediately inside the draft-tube exit. The 
shapes and sizes of these pieces as well as the nature of the broken material 
indicated that they were from the original half cones. No sections of the }-in. 
steel plate from the surface of the original cones were found in these two tubes. 
The draft-tube exit of Unit 1 was clear because approximately 75% of the 
original cone concrete in this tube was still in place. In this tube it was an easy 
matter to chip out the remains of the old cone because the concrete was of 
inferior quality and most of the steel plate had washed away. The steel shells 
of the original cones were filled through a small opening at the top; this was 
done despite a considerable handicap caused by leakage water from the penstock 
and the,threat of flood-water conditions in the tailrace. 


Expanding grout 


20 ft 3.75 In. 


(b) Section A-A 


(a) Pyramid 
reinforcing 


Fie. 12.—D1acram Snowine Reinrorcina Bars 


The new pyramids were installed by segregating (with sandbags) and drying 
the part of the draft tube involved. The original concrete surface was thor- 
oughly cleaned of scum and was carefully chipped to clean, sound material 
after which No. 11 reinforcing rods were grouted at least 4 ft deep into the old 
concrete (Fig. 12). Wooden forms were lined with an absorptive form liner 
# in. thick in order to produce a smooth, hard, dense surface on the finished 
pyramid’ Forms were not lowered into the draft tube until a few hours before 
pouring began because of the wet and humid condition. After the forms were 
assembled, several large light bulbs were installed inside them in order to help 
keep the absorptive liner dry. Concrete for the pyramids was a seven-sacks- 
per-cubic-yard design with an air-entraining agent; slump of from 2 in. to 4 in. 
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was specified. The average 28-day strength of concrete test cylinders was 
3,800 lb per sq in. A man was inside the form during the pouring process to 
work the concrete along the face of the forms with a spade as a precaution 
against honeycombing. A vibrator was used during the pouring, but special 
care was taken to prevent it from contacting the anchor dowels as the concrete 
in the bottom of the form had reached its initial set before the pyramid was 
“topped out.’’ The forms were removed 36 hr after concrete was poured, and 
the draft tube was flooded a day later. Each pyramid was allowed to cure for 
14 days with the unit shut down. From 14 days to 21 days after pouring, the 
unit was restricted to generation at best gate in order to keep turbulence in the 
draft tube ata minimum. After 21 days, the unit was released for unrestricted 
operation. 

Excellent results were obtained in producing a smooth surface on the pyra- 
mid and in using expanding grout around the periphery. The pyramids for 
all three units were installed at a total cost of approximately $41,000. 


Sew 


12000 16000 20000 24000 28000 32000 36000 40000 44000 48000 52000 
Turbine output, in horsepower 


Turbine efficiency, in % 


Fie. 13.—Tereine Erriciency, Unit 1 


In June, 1953, a second Gibson test was performed on Unit 1 for the pur- 
pose of comparing the performance of this unit after the pyramid was installed 
in the draft tube with the 1927 acceptance test. This test disclosed that the 
maximum turbine output was increased from 50,500 hp to 51,330 hp. This 
represents a gain of approximately 605 kw in generating capacity, which con- 
formed to the expec ‘ted results. The efficiency at maximum capacity was in- 
creased from 79.5% to 82.2%. Best gate operation was increased from 36,400 
hp to 39,000 hp with a corresponding increase in maximum efficiency of from 
90.3% to 91.1%. Although the efficiency decreased approximately 0.7% for 
gate openings producing less than 36,000 hp, its average increase was more than 
1% for those producing more than 36,000 hp. Both tests were based on a net 
head of 145 ft. The increased efficiency is in the range of operation where it 
will be most beneficial. The ‘kink,’ or abnormal drop in efficiency that oc- 
curred near 40,000 hp, without the modification, was definitely smoothed out. 
This is shown by a comparison of the two Gibson-test curves (Fig. 13). 

A comparison of the operating characteristics of Units 1, 2, and 3 at Martin 
Dam results in the following conclusions: 


An increase in generating capacity and efficiency was obtained by the 
installation of the half pyramid in the heel of the draft tubes (Fig. 10(b)), thus 
again indicating the value of hydraulic model studies. 

2. The installation was economically justifiable because the cost was o.iset 
by the value of the increased energy and capacity in approximately one year. 


748-24 


a 
5, 
+ 
~ 
| 
. 
4 
4 


Cart E. KinpsvaTer® Ric#arp R. Ranpowps, Jr.,"° Members, ASCE. 


* Prof. of Civ. Eng., Georgia Inst. of Technology, Atlanta, Ga. 
1 Hydr. Engr., Southern Services, Inc., Birmingham, Ala. 


—Mr. Woodruff’s discussion is a necessary complement to the writers’ descrip- 
tion of the events leading to the latest attempt to remedy the undesirable 
characteristics of the original draft tubes in the Martin Dam power plant. The 
details of design and the methods used in constructing the remedial structure 
are of particular interest because they were devised in the light of the un- 
fortunate experiences associated with the modifications made in 1928. The 
disintegration of the original remedial structures, revealed when the tubes were 
unwatered, is sufficient evidence of the necessity for the meticulous care with 
which the latest construction was accomplished. 

Mr. Woodruff reports that concrete remnants of the original heel cones were 
of inferior quality. The poor quality of the concrete used was doubtlessly 
indirectly related to the difficulty of placing the backfill under the outer steel 
shell. The combination of an infirm backfill and an insufficient anchorage 
apparently accounted for the failure of the steel armor under the unsteady 
dynamic loads resulting from turbine discharge. The simplicity of the new 
pyramid-shaped structures greatly facilitated the construction described by 
Mr. Woodruff. It is believed, furthermore, that the omission of a steel armor 
and the concentration on the production of superior-quality concrete which he 
describes shows a proper recognition of the necessity for fitting materials and 
construction methods to the conditions of the job. 

The writers were particularly pleased that comparative efficiency tests were 
made on Unit 1 before and after the remedial construction was completed. 
Prototype verification of conclusions based on hydraulic model studies is 
admittedly viewed with considerable relief by the researcher. It is true that 
the basic principles of similitude are as indisputable as the most fundamental 
principles of engineering mechanics. Unfortunately, however, the require- 
ments for complete similitude are not always possible to fulfil in the laboratory. 
It follows that a considerable amount of judgment, based on experience as 
well as theory, is an important component of every model study. The omni- 
scient Leonardo is quoted as a final warning: ‘‘Experiment never errs; only 
your judgment errs in anticipating results which experience does not confirm.” 

In connection with Mr. Woodruff’s discussion of the prototype tests, it is 
pertinent to emphasize the impracticality of a direct comparison of the results 
of the turbine efficiency tests with the results of the hydraulic model tests 
reported in the paper. Sufficient evidence is believed to be available to sub- 
stantiate the use of “diffuser efficiency”’ as an indication of the relative worth 
of different draft tubes. The writers do not mean to imply, however, that the 
draft-tube component of the power efficiency of a turbine can be determined 
from Eq. 5. For this investigation specifically there is little reason to believe 
that general draft-tube characteristics need be determined from a scale model 
of the complete machine. Indeed, it was the simplicity of the laboratory 
apparatus used by the writers that was unique in comparison with previous 
studies, 

It was observed following the derivation of Eq. 5 that » is an adequate 
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measure of the energy efficiency of a diffuser only when the area ratio (A3/A;) 
approaches zero. It should also be acknowledged that the writers’ interpre- 
tation of Eq. 5 is inadequate as a means of determining the true energy efficiency 
of any diffuser unless the flow at the end sections is both uniform and axial. 
It follows that, because the effect of whirl was ignored in evaluating Eq. 5 
from the laboratory data, an error of indeterminate magnitude is inherent in 
certain of the computed efficiency values. This does not of, course, invalidate 
the procedure used to determine the relative merits of the various draft-tube 
modifications. 

In the derivation of Eq. 5, it was implied that the flow was both uniform and 
axial. For this condition, the average piezometric head at both Section 1 and 
Section 3 is represented by the piezometric head at the boundaries. This 
fact is important, because only at the boundaries are the piezometric heads 
readily measured. With whirl at Section 1, however, the piezometric head 
actually varies from a minimum at the center to a maximum at the wall. The 
magnitude of this variation depends on the magnitude and distribution of the 
velocity in the cross section. Unfortunately, neither the velocity pattern nor 
the piezometric-head distribution can be measured by ordinary means. It 
follows, therefore, that the use of wall piezometers to evaluate A, results in an 
increasing overevaluation as the whirl increases. At Section 3, of course, the 
residual whirl is usually so small that it has negligible influence on the normal, 
hydrostatic pressure variation in the section. 

The second complication resulting from whirl at Section 1 concerns the 
denominator of Eq. 5. Thus, when the flow contains a whirl component, V; is 
correctly defined as the vector sum of the axial velocity and the tangential 
velocity in the section. For a given discharge, therefore, the average kinetic 
energy at Section | increases with increasing whirl. The exact magnitude of 
this influence is also indeterminate. It was expedient and consistent with the 
objectives of the investigation, therefore, to evaluate V, as the average axial 
velocity, Q/ A. 
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Discussion of 
"TRANQUIL FLOW THROUGH OPEN-CHANNEL CONSTRIC TIONS" 


by Carl E, Kindsvater and Rolland W, Carter 
(Proc, Paper 467) 


Emmett M. Laursen," A.M. ASCE, and Artuur Tocu,” J.M. ASCE.— 


11 Research Engr., Iowa Inst. of Hydr. Research, State Univ. of Iowa, Iowa City, Iowa. 
12 Research Associate, Iowa Inst. of Hydr. Research, State Univ. of Iowa, lowa City, Iowa. 


Within the limitations of the range of geometries used in the model investiga- 
tions, the coefficients presented by the authors should aid greatly in the in- 
direct determination of streamflows. The scatter of + 5% in Fig. 10 indicates 
a higher degree of accuracy for the discharge coefficient than is likely to be 
obtained for Ah and m in the field. Because the limitations recognized by the 
originators of such solutions as this are too often misunderstood or ignored, two 
major limitations on the applicability of this method of determining stream 
discharges should be mentioned. 

If the stream banks between zones 2 and 3 are covered with brush and trees, 
the water-surface elevation on the downstream side of the constriction may 
not be equal to that in the constriction. Depending on the density of the 
vegetal screen, the elevation of the water surface in zone 3 rather may more 
nearly approximate that elevation in zone 4. Although laboratory experiments 
with screens of known solidity ratio could determine this effect, in practice the 
difficulty of determining the density of a brush and tree cover would still 
remain. The companion work" to the paper might be used to obtain Ah in 


3 Backwater Effects of Open Channel Constrictiohs,"" by Hubert J. Tracy and Rolland W. Carter, 
Transactions ASCE, Vol. 120, 1955 


this case, but in using the curves therein, it is suggested that m be determined 
from the downstream discharge ratio. Because the slopes of the main channel 
and the flood plain would not be the same, the conveyance ratio should not be 
used. 

If the bed of the stream is erodible, scour in the constriction may occur 
during a flood and the scour hole may be completely or partly refilled by the 
time the elevation of the bed of the stream is determined. Considerable error 
could thus, of course, be caused in the estimation of the depth of flow. An 
approximate indication of the order of magnitude of this effect has been made" 


4 Discussion by Emmett M. Laursen and Arthur Toch of ‘‘River-Bed Scour During Floods,"’ by E. W. 
Lane and W. M. Borland, ibid., Vol. 119, 1954, p. 1084. 


by the writers. 

The writers wish to express their admiration for the careful experimental 
work which forms the basis of this paper. It is agreed that the analysis, 
admittedly empirical, is the only one feasible; much more information is 
needed about the fundamentals of turbulent flow before a more rational 
analysis is practicable. 

Frep W. BraispeE.t,'* M. ASCE.—When the results of a comprehensive 


18 Project Supervisor, U _8. Dept. of Agriculture, Agri. Research Service, St. Anthony Falls Hydr. 
Lab., Univ. of Minnesota, Minneapolis, Minn. 


series of generalized tests are summarized in the few curves presented by 
Messrs. Kindsvater and Carter, it is difficult to realize the tremendous effort 
and large number of tests that are required to obtain these results. The writer 
and Charles A. Donnelly completed (in 1950) a generalized study of the 
eapacity of the box inlet drop spillway ;'* the methods of analysis and presenta- 
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16 “Hydraulic Design of the Box Inlet Drop Spillway,” by Fred W. Blaisdell and Charles A. Donnelly, 
paeete Paper No. 8-B, St. Anthony Falls Hydr. Lab., Univ. of Minnesota, Minneapolis, Minn., January, 
tion closely paralleled those used by the authors. Because of this experience, 
the writer can fully appreciate what is meant by the statement (under the 
heading, “Simple Constrictions in Rectangular Open Channels: Experimental 
Verification of the Analysis’) that ‘After much study, a procedure was 
adopted ***.”’ The benefit is well worth the effort for, in presenting their 
results, the authors have done much of the work required for the practical 
application of the results of their research. 

As noted by Messrs. Kindsvater and Carter, the discovery that variations 
in the geometrical shape of the channel cross section ean be expressed in terms 
of conveyances is certainly an interesting and valuable contribution. It is 
to be regretted that values of m less than 0.20 were not included in the test 
schedule. It would have involved only a little more work and would have 
eliminated questions that always arise when curves are extrapolated beyond 
the range of the tests. It may be that these low values of m are assumed to be 
outside the range of practical application. However, there may arise practical 
situations in which values of m will be less than 0.20. It has been the writer’s 
practice to extend the range of laboratory tests well beyond the “practical” 
limits set by those experienced in the actual design of structures. Even so, 
subsequent practical applications outside the extended range of test variables 
have been frequent. 

Tests at values of m less than 0.20 would have been particularly helpful in 
defining the effects of angularity, corner rounding, and guide walls. In the 
case of the adjustment curves for angularity, a large part of the curvat ire is at 
values of m less than 0.20. Although the curves are shown to be solid in Fig. 8, 
One must assume (based on the table in Fig. 10) that they must huve been 
drawn without the benefit of experimentally determined points. 

The comparison of the adjusted and observed values of C presented in 
Fig. 10 provides a crucial test of the analytical methods. The agreement 
shown is excellent. Although it is undoubtedly possible to obtain less spread 
in the data for individual curves, greater spread must be expected as adjust- 
ments are added. Overcompensating benefits include an expanded range of 
application and an increased reliability resulting from the greater volume of 
data required to define the coefficient and adjustment curves. 


Pin-Nam Lin.'’—The discharge coefficient has been virtually detined by 

7 Asst. Prof., Dept. of Civ. Eng., Colorado Agri. & Mech. College, Fort Collins, Colo. 
Messrs. Kindsvater and Carter in two ways. It is defined in one case by the 
first dimensionless group in Eq. 6b and in the other by Eq. 10. Because hy 
and a, V*,/2 g in Eq. 10 have opposite signs it is possible that in many cases 
the two values of C are approximately the same. In general, the two defini- 
tions will yield different values of C; the case of channels having considerable 
tree growth and brush, cited by the authors, is an example. The apparent 
discrepancy, however, appears to be the result of an error in presentation 
because, if the first definition were to be considered, Eq. 6a should have in- 
cluded the coefficient a, as well as a variable to indicate the roughness of the 
channel. For a given boundary geometry at a given locality, Ah will vary not 
only with Q and y2 but also with the channel roughness and the velocity 
distribution upstream of the constriction. When C is defined by Eq. 10, it is 
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reasonable that C will be expressed by Eq. 12 or Eq. 17. 

Flow through channel constrictions at high values of m is accompanied by 
considerable surface curvature. As a result of the surface curvature, the 
accéleration in a vertical plane may be considerable. Analysis of a flow 
that is no. longer “gradually varied”’ was first conducted by Boussinesq.'*% 


18 Steady Flow in Open Channels: The Problem of Boussinesq,”’ by C. Jaeger, Journal, Inst. of C. E., 
London, England, February, 1948, pp. 339-349. 


'9*Fluid Mechanics for Hydraulic Engineers,"’ by H. Rouse, McGraw-Hill Book Co., Inc., New York, 
N. Y., 1938, pp. 301-307. 


Because the pressure distribution is no longer hydrostatic in the section where 
the surface curvature is pronounced, coefficients should also be introduced to 
account for the effect of pressure distribution on the static-head terms in the 
energy and momentum equations. Because the surface curvature at section 
2 may be considerable, it would be more appropriate to write Eq. 7 as 


V2, 


In tranquil flow, the static head above the bottom of the channel is usually the 
major part of the specific head so that—depending on the surface curvature— 
8 may be an important coefficient. In fact, because the flow at section 2 is 
markedly convergent the velocity distribution in this section should be ap- 
proximately uniform, and thus a2 probably does not differ greatly from unity; 
that is, a, is often of secondary importance. From Eq. 21, one would be in- 
clined to include the coefficient 8 in Eq. 9. It is believed that this step is 
probably not necessary in the present case, in which the main objective is to 
compute the discharge on the basis of experimental coefficients of discharge. 
One may consider Eq. 10 or Eq. 13 as a formula defining C, which is expressed 
by Eq. 12 or Eq. 17; the effect of surface curvature is included in Eqs. 12 and 
17. 


in which 


If, however, Eqs. 9 and 10 are regarded as rational deductions from the 
equations of energy and continuity (that is, Eqs. 7 and 8), the coefficient 6 
should be included in the energy equation as indicated by Eq. 21. Boussinesq 
assumed” that the curvature in a vertical section varied linearly from the 

Technische Hydraulik,’’ by C. Jaeger, Birkhauser, Basel, Switzerland, 1949, pp. 116-122. 
surface to the bed and obtained, for the case of uniform velocity distribution in 
a vertical direction, the formula: 


Assuming further that the velocity distribution is uniform in section 2, there 


results 
J, (2 +2). 
, 


748-29 


> 


i | 
1 
| 4 
| 
> 
i 
B 
| 
| 
| 
| 


(v 
2 


In general, C. Fawer® showed that, by assuming the curvature at any level 


a “Etude de quelques écoulements permanents & filets courbes,” by C. Fawer, La Concorde, Lausanne, 
Switzerland, 1937. 


as given by (z/y)" d*y/dz*, Eq. 24 could be extended to yield 


n+1lgy 


Naturally, both Eqs. 24 and 25 are far from being exact. These expressions 
are, nevertheless, at least qualitatively correct. Because the surface at 
section 2 is concave d*y/dz? is positive, and therefore 8 > 1. This result may 
be anticipated if it is noted that the centrifugal effect of flow at section 2 is to 
increase the local pressure beyond that predicted by the hydrostatic pressure 
distribution. For a given value of m, different cases of flow through a given 
type of constriction should be geometrically similar. Asa result, corresponding 
to a given value of m, the dimensionless parameter y(d?y/dz*) should be ap- 
proximately constant. On the basis of Eqs. 24 and 25, several qualitative 
observations can thus be made concerning 8 for flow through a given type of 
constriction: (1) Even if m is kept constant, 6 still increases with the Froude 
number. (2) Regardless of the magnitude of m, 8 approaches unity as the 
Froude number approaches zero. (3) The term @ increases with m. These 
observations indicate clearly that the magnitude of 6 depends not only on the 
surface curvature of flow but also on the Froude number. 

In order to gain some idea regarding the influence of the coefficient 8, the 
possible effect of 8 on the ratio (y: — ys)/Ah will be cited briefly. Because 
the value of 6 at section 2 is greater than unity, the total energy per unit 
volume of water at section 2 is greater than that indicated by (y2 + az V*,/2 g). 
The value of Ah computed by the energy equation without taking 8 into 
account will therefore be less than the actual value of Ah. The effect of 6 is 
thus to increase the value of Ah. However, the quantity (y: — y,), being 
governed essentially by conditions in sections 1 and 4, may be considered 
independent of 8. Therefore, the resultant effect of 8 is to reduce the ratio 
(y: — ys)/OA. In the light of the three observations previously made regard- 
ing 8, the following remarks concerning (y: — ys)/Ah can be made: 


a. For a given value of m, (y: — ys)/MhA decreases as the Froude number 
increases. 

b. For a given value of the Froude number, the decrease in (y; — ys)/AA 
is greater as m increases. 


Both of these trends appear to be supported by experimental data.” Con- 


Backwater Effects of Open Channel Constriction,” by Hubert J. Tracy and Rolland W. Carter, 
Transactions, ASCE, Vol. 120, 1955, igs. 6-8. 


sequently, within the range of data presented by Messrs. Kindsvater and 
Carter, it seems certain that a more rational analysis of the problem should 
include a consideration of £. 
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Cart F. Izzarp,* A.M. ASCE.—The placing of a constriction in an open 
% Chf., Hydr. Research Branch, Bureau of Public Roads, U. 8. Dept. of Commerce, Washington, D. C- 
channel inevitably results in an increase in the upstream stage over the stage 
which existed for the same discharge in the unconstricted channel. In Fig. 11 


Fig. 11.—Derrtnition Sketcu oF BACKWATER 


there is illustrated the change in profile from flow at a uniform depth y, to a 
typical backwater flow. Also shown are the lines of total head obtained by 
plotting the velocity head above the free water surface, assuming that one can 
estimate the kinetic energy in the flow. The depth y, is the increased depth 
which Messrs. Kindsvater and Carter have measured at a section a distance b 
upstream from the upper edge of a constriction whose width normal to the 
direction of flow is b. To generalize the analysis, the channel is shown to have 
a slope So which actually was zero in the authors’ experiments. 

Using the nomenclature adopted by the authors, section 3 is taken at the 
downstream edge of the constriction. The drop in the total-head line from 
section 1 to section 3 can be considered to be equal to the normal resistance 
loss plus the loss h. caused by the contraction of flow and is probably small 
compared to the total loss. This flow, however, has to expand and in so 
doing additional energy is released that is represented by the difference h, 
between the total head at section 3 and the total head with uniform flow. 
The sum 


thus represents the total energy loss in excess of the normal friction loss between 
section | and a section downstream where the flow is again at normal depth. 
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From Fig. 11 it is apparent that A, is also equal to the change in total head at 
section 1. 
The energy equation between section 1 and section 3 then becomes 
29 29 
Where a cross section is irregular, a suitable value for the weighted velocity 
head should be substituted. Because h, results from a contraction of flow 
one can assume that it is a function of the velocity head at section 3; the 
experiments, however, have shown that the water surface is not level trans- 
versely at section 3. To avoid this complication one can assume that the area 
at section 3 is taken at the normal depth, relying on the head-loss coefficient 
determined empirically to correct for the difference. Using the symbol Vas 
to designate the velocity at normal depth in section 3, 


Substituting this value into Eq. 27 and transposing terms, 


V2, V2, 


+ Ky 


29 


From the continuity equation, Eq. 29 can be expressed as 


from which 


If A,/A,; approaches unity as it will for relatively small quantities of back- 
water, an approximate equation for Ky is 


Backwater can then be estimated as 


Yi — yn = Ky 29 
provided values of Ky, for a given type of constriction are known. After a 
trial value of y, is determined from Eq. 33, a more accurate value can then be 
obtained using Eq. 30 because a close approximation of A, will then be known. 

Experiments have been described'* from which it has been possible to test 
the applicability of Eqs. 29 to 33. In those tests steady flow was first estab- 
lished with no constriction, and an equivalent of normal depth was recorded. 
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The flume used had a level floor, but the change in depth between section 1 
and section 3 was usually very small in relation to the change in depth between 
the same sections with the constriction in place. The most extensive data 
were for a vertical constriction with a symmetrically placed rectangular opening. 
In Fig. 12 it is shown that K, varies systematically with the contraction ratio 
m and to a lesser extent with the parameter V*,3/(2g yn). In Fig. 12 only 


16 


Value of K, 
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Fig. 12.—Backwater Coerricient as A Frnerion or Vetocity Heap Revative 
To Norma. Deptn (n = 0.012, L = 0.02 ft) 


tests made with L = 0.02 ft and n = 0.012 are shown. Thecorrelation is good 
except for a few tests, in which the surface profile for flow with no constriction 
appeared to be irregular. A more significant plot of the data is shown’ in 
Fig. 13(a) which is derived from the curves drawn in Fig. 12. From Fig. 13(a) 
it is indicated clearly that K, is primarily a function of m, the variation with 
the velocity-head parameter being of minor importance; similar and slightly 
lower curves could be drawn for other values of L. Roughness of channel ap- 
péars to affect K, slightly, but this difference may be the result of the increased 
head required to obtain flow on a level floor. 

A few runs were made with other types of constrictions. The computed 
values of K,, for two types are shown in Fig. 13(b) in comparison with the curve 
for the rectangular opening. The 45° wing wall had a constant dimension 
w = 0.25 ft so that w/b was variable with the value of b and therefore with m; 
this accounts for the tendency of the curves for the wing-wall type to approach 
the curve for a rectangular opening when w is relatively small, and to approach 
the curves for the rounded embankment for the case in which w is relatively 
large. The radius at the toe of the | on | rounded embankment was not 
constant, this radius varying from one third to one half of the distance across 
the constriction between the toe of the slope on each side. 
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Fig. 13.—Backwater Coerricient as a FuNction oF Contraction Ratio (n = 0.012) 


These analyses, based on meager experimental data, strongly suggest the 
possibility—using a tilting flume—of determining values of K, that would be 
applicable to standard types of bridges and that would be independent of 
channel roughness. If so, backwater could be computed by using Eq. 33. 
By applying this method the contracted depth y; at the downstream side of 
the abutment would not be given, but this value is not important to the bridge 
designer. Actually y, — ys is-a function of V?,3/2 g similar to Eq. 33. 

The foregoing discussion has been directed at the problem of estimating 
backwater ¢aused by constricting a channel carrying a known flow at normal 
depth, which is the problem of the highway engineer. Messrs. Kindsvater and 
Carter, however, were confronted with the reverse problem of estimating 
discharge from measurement of water levels in the vicinity of a channel con- 
striction. Because the depth y, is sensitive to small changes in geometry, it 
Was necessary to devise a system which would reflect these changes; this they 
have done with remarkable ingenuity. 

Performing the tests in a channel with a level floor greatly simplified the 
experimental procedure but left some question as to the applicability of the 
discharge coefficient to flow in a sloping channel, particularly when h, is large 
in comparison to Ah. The basic equation defining C (Eq. 13) includes a term 
for head loss in the approach reach computed by the geometric mean of the 
energy slopes at sections | and 3 and the energy slope through the contracted 
reach. The evaluation of C would have been simpler if h; had been computed 
solely for the friction loss in the unobstructed channel over the length b + L. 
All the energy losses attributable to the geometry of the constriction would 
therefore be included in the discharge coefficient. Thus, one would avoid the 
question as to whether or not the friction loss is correctly evaluated by the 
geometric mean of the energy slopes in the case in which the transition is not 
gradual] but abrupt. 

For computing discharge for an actual bridge, the normal friction loss could 
be checked against the energy slope in an unobstructed reach of the same cross 
section and should be approximately equal to the average bed slope. Thus, 
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the uncertainties involved in estimating Manning’s n—required to compute the 
friction loss—would be avoided. It would still be necessary to use n in es- 
timating the contraction ratio but, because n appears in both numerator and 
denominator, the effect of error would be minimized. Consequently, the 
evaluation of C in Eq. 13, using Ay as the normal friction loss, would simplify 
computation in the field. Further simplification could be attained by defining 
m as K,/Kg instead of 1 — K,/Kg. The obvious term for m would then be 
“conveyance ratio.” 

The problem of estimating discharge from field measurements requires 
that the elevations h, and h; must be measured at certain points to determine 
Ah. In estimating backwater y; — yx, it is not necessary to compute ys, if 
K, in Eq. 33 is known from experimental data. The computation is simple 
and gives results which appear to be as consistent, if not more so, than those 
obtained by computing Ah from Eq. 13 and then applying corrections to get 
Yi — Yn. 

In a discussion™ the writer demonstrated that y:/y, was a function of 


™ Discussion by Carl F. Izzard of ‘Backwater Effects of Open Channel Constrictions,"’ by Hubert J. 
Tracy and Rolland W. Carter, Transactions, ASCE, Vol. 120, 1955 


V?,.3/(2g 9.) and the parameter m for a given type of constriction. By sub- 
tracting unity from the ordinate scale and multiplying both scales by yp, 
these quantities take the form of Eq. 33 in which K, becomes the slope of a 
straight line approximating the curved lines. The reason for the consistency 
of the data in relation to the velocity-head parameter based on normal depth 
in the contracted opening is that the total head loss caused by the contraction 
h, (as defined in Fig. 11) is referred to total energy at normal depth. 

The use of Eq. 33 together with charts similar to Fig. 13(6) is believed 
preferable to the use of another chart?® because Eq. 33 gives the increment in 

% Ibid., Fig. 9. 
depth and not a depth relative to a normal depth which might prove difficult 
to define for irregular channels. More experimental data are needed, prefer- 
ably from a sloping flume, to define K, for standard types of bridges 

Cues.tey J. M. ASCE.—The work of Messrs. Kindsvater and 


c 26 Prof. of Hydraulics and Structural Eng. and Head, Dept. of Civ. Eng., State Univ. of Iowa, lowa 
ity, Iowa. 


Carter should do much to clarify the thinking of engineers on a problem that 
has often proved puzzling. The tests reported not only provide useful in- 
formation but also point the way for future experimentation that will be even 
more valuable. 

The authors properly emphasize the importance of accounting for the con- 
traction of the jet. They omitted one factor, however, that can have a con- 
siderable effect on the contraction—friction along the upstream face of the 
obstruction ; the corresponding orifice effect is recorded in standard text books 
on hydraulics.2”7 Messrs. Kindsvater and Carter should have included some 


oes. 1 apamedl by H. W. King and C. O. Wisler, John Wiley & Sons, Inc., New York, N. Y., 3d Ed., 


tests with variable roughness in this location inasmuch as they did vary the 
bed roughnesses. 

In their future tests, the authors should give special attention to the con- 
traction to be expected under various conditions. The right of way along 
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the fill section approaching a bridge is apt to be kept cleared whereas the 
flood-plain bottom beyond the right-of-way lines is often overgrown with un- 
derbrush and trees. Hence, the flow parallel to the fill approaching the bridge 
will attain a higher velocity than it would under the flow-net assumption of 
equal (or nonexistent) friction throughout the approach area. Entering the 
main flow transversely just above the bridge opening, this current of high ve- 
locity will cause more contraction than would otherwise exist. Not only is 
this condition almost universal but it is often accentuated by the presence of 
borrow pits or large ditches along the side of the fill. 

The designation of some formulas as “rational” and others as ‘‘empir- 
ical” seems unfortunate. Because all formulas representing physical phenom- 
ena are empirical, what is probably meant is that some have more claim to 
rationality than others.2* The rationality of the Weisbach formula is not at 


**’The Science of Mechanics,” by Frnst Mach, Open Court Publishing Co., Chicago, IIL, 5th Ed., 
1942. 


all appealing but, because it includes two terms with different exponents (and 
hence is similar to a power series), it can be fitted more closely to the experi- 
mental data than can other formulas that have a more logical basis but only 
one term.’ A close examination of the Nagler formula’ shows that it has at 
least as good a claim to “rationality” as has, for example, the application of 
the flow net to open-channel problems. The authors’ introduction of the con- 
cept of contraction in the study of open-channel flow is not new*” but, when 

» An Experimental Study of the Flow of Water through Transitions in Rectangular Open Channels,” 
by George B. Lyon, thesis presented to the State University of lowa, lowa City, lowa, in February, 1942, 
in partial fulfilment of the requirements for the degree of Master of Science. 
correlated with systematic experimental studies, it seems to bring order out 
of chaos. One wonders whether Messrs. Kindsvater and Carter could predict 
the experimental results obtained® by Mr. Lane by ana priorievaluation. Per- 
haps their suggestion that the degree of jet contraction should be proportional 
to the ratio m can be provided with “‘adjustments” to take into account the 
existence and relative elevations of side borrow pits, ditches, or cleared right 
of ways. 

Cart E. Kinpsvarer,” M. ASCE, anp W. Carrer,” A. M. 


* Prof. of Civ. Eng., Georgia Inst. of Technology, Atlanta, Ga. 
* Hydr. kEngr., Geological Survey, U. 8. Dept. of the Interior, Atlanta, Ga. 


ASCE.—The subject of this paper was the first phase of a comprehensive in- 
vestigation of the flow of water through bridge openings. Since the delivery of 
the paper at New Orleans, La., in 1952, two additional reports on related 
research have been published. In the first of these” the study of discharge 


#2” Computation of Peak Discharge at Contractions,"’ by Carl FE. Kindsvater, Rolland W. Carter, and 
Hubert J. Tracy, Circular Nv. 284, Geological Survey, U.S. Dept. of the Interior, Washington, D. C., 1953. 


characteristics is extended to include additional boundary conditions; the 
second" deals with the flow pattern downstream from the constriction. As the 
basis for these and other investigations still in progress, the writers’ paper is 
appropriately complemented by the excellent discussions submitted by Messrs. 
Laursen and Toch, Blaisdell, Lin, Izzard, and Posey. 

Messrs. Laursen and Toch refer to certain limitations on the applicability 
of the writers’ solution. They observe correctly that the relationship between 
the level of the water in zone 3 and that in zone 2 depends somewhat on the 
density and distribution of the vegetal cover on the banks of the stream. It 
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should be emphasized, however, that the experimental coefficients are based on 
measurements made at zone 3A. Jn terms of natural channel conditions, this 
zone is on the downstream face near the waterway end of theembankment. It 
is unlikely that the influence of trees and brush between this point and zone 2 
would ordinarily be significant. a3 

Mr. Lin’s remarks concerning contradictory definitions of C appear to be 
based on a misunderstanding of the purpose of the dimensional analvsis leading 
to Eq. 6b. In the section entitled “Simple Constrictions in Rectangular Open 
Channels: Selection of the Independent Variables” the resistance characteris- 
tics of the upstream channel are recognized, but they are omitted from the 
dimensional analysis in anticipation of their subsequent treatment in the dis- 
charge equation. Thus, the second dimensionless group in Eq. 66 is described 
as a ratio that is “proportional to a coefficient of discharge.” It is, in fact, 
simply an indication of the existence of a general parameter having the form of 


a Euler number. It was the purpose of the section following Kq. 6 to formulate | 
the coefficient of discharge expressed by Eqs. 10 and 12. A brief justification of $ 
the writers’ method of accounting for h; and a, V?,/2 g is contained in the para- . | 
graph following Eq. 9. | 
Mr. Lin describes some of the complications that result from the application ces 
of the one-dimensional energy and momentum equations to a section in which . | 
the pressure is nonhydrostatically distributed. It is usually possible to avoid | 


these complications by applying the equations to zones of parallel flow; that is, 
to sections in which the normal acceleration—and, therefore, the normal piezo- : | 
metric-head gradient—is zero. This was accomplished in the writers’ study by a | 
defining section 2 quite arbitrarily as the vena contracta. Actually, of course, er | 
this qualification is largely academic. Regardless of where section 2 is located, A ‘| 
both a» and B» ordinarily will be functions of the independent variables upon . 
which the experimentally determined values of C depend. It follows that the 
influence of a nonhydrostatic pressure distribution, in addition to the influence if 
of a nonuniform velocity distribution in the constricted section, is contained in : 
the coefficient of discharge defined by Eq. 13. 
The writers are grateful for Mr. Biasidell’s recognition of the large amount iP +| 
of work and the many tests required to establish empirical solutions for complex 
fluid phenomena. Problems involving open channels, particularly natural 
channels, typically depend on many independent variables. It is seldom pos- 
sible to consider all the variables in a laboratory investigation. The decision to 
limit the experiments to values of m greater than 0.20 was one of the conse- : 
quences of the necessity to restrict the scope of the study to the conditions most a 
frequently encountered in practice. Nevertheless, it should be recognized that ; 
the tests were conducted at a relatively small scale and therefore that values Ee 
of m less than 0.20 would result in very small changes in water-surface level and 
quantities such as Ah would be subject to large experimental errors. The | 
writers agree with Mr. Blaisdell that, from the viewpoint of the researclier, so- ; : 
called “practical” limits are often frustrating. 
Much of Mr. Izzard’s discussion is concerned with the determination of the q 
backwater due to an open-channel constriction. This is one of the objectives 4 : 
of the comprehensive investigation in which the writers are engaged. An 
analysis and much additional experimental data pertaining specifically to that | 
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problem are therefore contained in a subsequent paper." 

Mr. Izzard is concerned about the fact that, in all the writers’ tests, the 
bottom of the experimental flume was horizontal. The significance of this 
limitation is that a small part of the change in depth in a level flume is due to 
the acceleration which results from the unbalanced boundary resistance force. 
For the conditions of this investigation the influence of bottom slope is believed 
to be insignificant. Mr. Izzard also proposes a simplification in the method of 
computing h;. In his words, “All the energy losses attributable to the geome- 
try of the constriction would therefore be included in the discharge coefficient.” 
Actually, he would be adding one variable, the conveyance of the constricted 
waterway, to those already involved in the evaluation of the discharge coeffi- 
cient. That this variable has been effectively removed in the writers’ solution 
is demonstrated by the lack of correlation between C and n throughout the full 
range of tests reported. In connection with both of Mr. Izzard’s criticisms, it 
should be observed that the contracted-opening method of determining the 
discharge in natural channels is not feasible when the boundary resistance loss, 
hy, is « large part of Ah. 

Mr. Posey deplores the omission of certain natural channel variables, in- 
cluding constriction roughness as well as borrow pit, drainage ditch, and right- 
of-way characteristics. In a subsequent report many variables were con- 
sidered in addition to those involved in the present study. These variables, 
all concerned with natural bridge waterways, resulted from sloping embank- 
ments, a variety of abutment forms, piles and piers, and partial bridge sub- 
mergence. Thus, of the “almost unlimited” number of variables acknowledged 
in the “Introduction,” a large number have been considered. It is possible that 
experience in the application of the proposed computation procedure will in- 
dicate the necessity for experiments involving additional variables. In the 
meantime, growing evidence of the successful application of the writers’ data 
to the computation of flood flows in natural waterways is encouraging. 


748-38 


tag 
| 
| 
‘| 


| 


Discussion of 
“THE PRESENT STATUS OF RESEARCH : 
ON SEDIMENT TRANSPORT” 
by Ning Chien po) 
(Proc. Paper 565) { 
» 


HSIN-KUAN LIU.!—To make a brief summary of what has been studied and 
found in the field of sediment transport is not an easy task, because numerous Pe 3 
papers dealing with the mechanics of this subject have been published. It is a 
even more difficult to bring conclusions from various studies because of the 
widely spread opinions among them. The author certainly has done a timely 
summary which helps the advancement of research in the field of sediment 
transport. 

In discussing the initiation of ripples, the author introduced the Writer’s 
result as shown in Fig. 2 Curve II (Fig. 7 in this discussion). The result is 
based upon the assumption that ripples are due to instability of the interface 
(or boundary) between the flow and the movable bed. The author argued that 
since the sediment rate is not large enough to create any difference in density . 
between the clear water and the sediment-laden water, ripples cannot be as- ; i 
| 


sociated with the instability at the interface of two flows of different densities. 
The writer suggests the following as a possible explanation of the mechanics 
of ripple formation. 

Due to difference in velocity between the flow and the movable bed, there 
exists an interface having a high velocity gradient. Furthermore, such an 
interface is very unstable because of this high-velocity gradient. Since the 
bed is movable, the bed configuration is readily affected by temporal varia- 
tions of flow pattern, if the local flow velocity is large enough to transport the 
bed material. The next question is what causes this temporal variation of 
flow pattern which results in the instability of the interface. This question 
can be answered, at least in part, by the analysis which Rouse used in explain- 
ing the instability of flow at a hypothetical surface of discontinuity. (94) At 
any instant the streamline at the interface may deviate from a straight line 
due to a disturbance which likely exists in the flow. The effect of the disturb- 
ance upon the flow pattern is dependent upon the instability of the interface 
(or boundary). According to the flow net, the local velocity will decrease as 
the neighboring streamlines diverge, and increase as they converge. If the 
local velocity of the flow is large enough tc transport the sediment, the peri- 
odic change in velocity along the interface will result in a change of bed con- 
figuration. Exner (22) has established a differential equation for sediment 
transport in two-dimensional flow. That is Eq. 40 indicates that the variation 


=0 (40) 


1. Asst. Prof., Dept. of Civ. Eng., Colorado A & M College, Fort Collins, 
Colo. 
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where, 1 is the bed elevation, 
t is time, 
K is a factor relating sediment discharge to flow velocity, 
v is the flow velocity near the bed and is a function of distance 
along the flow direction, 
x is the distance in the downstream direction. 


‘ of bed elevation is due to velocity variation with respect to distance in the 7 
flow direction, 
Es An index indicating the instability of the interface can be expressed by a 
3 dimensionless parameter, py , which is a special form of Reynolds number. 
On the one hand the greater the shear velocity, which means a higher velocity 
gradient at the interface, the more pronounced will be the instability of the * 
~ a interface; and the larger the sediment size, the more pronounced will be the 
: tendency to produce eddies. On the cther hand, the higher the viscosity, the 
ae more the internal shear will tend to oppose further differences in velocity. 
Be The capability of the flow to transport the sediment and change the bed con- 
figuration is dependent upon the scouring force of the flow and the resistance 
force of the sediment. The scouring force can be expressed as (24) 


2 
Fi=Cig Wd*(apVy") (41) 
the resistance force of the sediment to the flow is (24) 
2 
Fe= (42) 


where Fy is the erosive force of the flow, 
Fa is the resistance force of the sediment, 
D is the sediment size, 
P¢ is the fluid density, 
V« is the shear velocity of the flow, 
Vg is the fall velocity of the sediment, 
vis the dynamic viscosity of the fluid, 


C;,C2 are functions of pee and shape factor of the sediment. 


: Laboratory experiments indicate that ripples occur shortly after the begin- 
zi. ning of bed-load movement. Therefore, it can be assumed without serious 
error that, at the moment ripples occur, the scouring force is equal to the 
critical scouring force. 


F,= Fp (43) 


so that 


DV 
= function shape factor) (44) 


oa If shape factor is considered to be of secondary importance, then Eq. 44 can 
d be written as 4 
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_ At the moment ripples occur, two conditions must be satisfied: (a) the 
flow must be strong enough to transport the sediment; (b) the interface be- 
tween the flow and the movable bed must become unstable. It should be noted 
that oye as mentioned above is a parameter which serves two functions: 

(a) it indicates the instability of the interface; (b) it correlates the flow veloc- 
ity at the grain level to the shear velocity in accordance with Karman-Prandtl 
theory, thereby defining the scouring force. By plotting data which are the 
average between smooth bed and first appearance of ripples in terms of 


Vs DV+« 


vy, versus = — Fig. 7 was obtained. Conditions (a) and (b) are both satisfied 
s 


when data pertaining to the flow and the sediment fall above the curve. The 
scatter of the data is due to such causes as error in measurements, difficul- 
ties of defining the first appearance of ripples, shape factor, and effect of 
mixture. 

The author pointed out that there are many bed-load formulas today. How- 
ever, their application to engineering problem is limited because of various 
reasons. The ones such as Einstein’s, Kalinske’, and Meyer-Peter’s are 
sound in theory, yet too complicated for many hydraulic engineers to under- 
stand. Furthermore, since there is no way to measure the actual bed-load in 
a large river, their applicability to a large river is uncertain. Many investi- 
gators are still searching for a better bed-load formula and also method of 
measuring bed-load. 

The writer’s experiment on sand wave movement, as will be explained 
below, seems to indicate the possibility of obtaining a type of bed-load formula 
in terms of sand wave motion, and of measuring the amount of bed-load trans- 
ported by investigating the sand wave motion. 

In performing the experiments the writer used a recirculating, open-chan- 
nel, laboratory flume 40 ft long by 1 ft wide by 2 ft deep. There were two 5-ft 
glass sections on one side and one 5-ft section on the other to facilitate obser- 
vation. Water was supplied from a floor-level sump of 8 ft by 2.5 ft by 2.5 ft 
through a centrifugal pump to a 4-in. feed pipe. Flow was measured through 
an orifice meter and returned to the supply sump over an adjustable tailgate 
at the downstream end of the flume. Sand was supplied at the upstream end 
of the flume by means of an elevator. The elevator was operated by a variable 
speed motor, chain, and sprocket-wheel connected to a screw-lift mechanism. 
Ottawa 20-30 sand was used as sediment. Sand transported to the downstream 
end of the flume was trapped in a bucket and weighed to verify the rate of sed- 
iment feed. The experiments were done in the St. Anthony Falls Hydraulic 
Laboratory of the University of Minnesota. The whole work was part of the 
writer’s doctoral dissertation under the guidance of Professor Lorenz G. 
Straub. 

The purpose of the experiments was to determine the quantitative relation- 
ship between sand wave motion and bed-load transport. It was found that ap- 
proximately 100 per cent of the bed-load is transported in the form of sand 
waves, except when there is an appreciable amount of saltation load in the 
stream. In these experiments, the discharge was kept constant at 0.4 cfs. The 
rate of sediment feed ranged from 600, 300, 150, 75, to 38.3 lbs per hour. The 
bed slope was varied from 0.0115 to 0.002. Mean velocity of the flow varied 
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number ranged from 0.97 to 0.512. 


(d = 0.690 mm.) 


c = 0.0008 


velocity. 


> as siderable amount of it was moving as saltation load. 


Results previously shown suggest that a formula for bed-load transport in 
a terms of sand wave motion, sediment properties, and flow properties might 
i possibly be obtained by extensive experiments. The formula, if obtained in 
A the future, could be used to verify existing bed-load formulas. Field meas- 
q urements of rate of bed-load transport could be made by observing the sand 
wave velocity during flood stage, and the bed profile could be taken by echo 


sounding devices. 


REFERENCE 


Inc., New York, 1947. 
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from 2.3 to 1.5 fps. Depth of flow varied from 0.174 ft to 0.267 ft. Froude 


The velocity of the sand waves was measured by two methods; (a) counting 
the number of waves passing a certain station per unit time and multiplying 
this by the wave length; (b) measuring the time required for a certain wave to 
travel a certain distance. The wave length was obtained by averaging the total 
length of several waves. It was found, as shown in Fig. 8, that the velocity of 
the sand waves for the case where q = 0.40 cfs with Ottawa 20-30 sand 


where c is the sand wave velocity and V is the mean velocity of flow. Inas- 
much as the wave velocity varies with a large exponential power of the mean 
velocity of the flow, the wave velocity is highly sensitive to change of the flow 


In studying the relationship between sand wave motion and rate of bed-load 
transport, the bed profile was measured by a point gage with an electric de- 
vice attached. As the point touched the wetted bed surface, a red light ap- 
peared on the gage. The depth of the flow was the vertical distance between 
water surface and the mean bed elevation. In order to compute the bulk vol- 
ume of the sand wave during each run, five bed profiles were taken along the 
flume. The lower boundary of the sand waves was considered to be the line 
joining the lowest points of the troughs. The average volume per sand wave 
was that of the total volume of sand waves above the lower boundary divided 
by the number of waves. The volume of sand waves per unit width was meas- 
ured from the bed profiles by a planimeter. Profiles of Run 2-2 are shown in 
Fig. 9. Results of the experiments are listed in Table I. The last two lines 
in Table I show both the actual and the computed amount of sediment trans- 
ported. It indicates that bed-load is entirely in wave forms. Translation of 
sediment along the bed is negligible. The computed amount of sediment in 
Run 2-1 is considerably below the actual amount of sediment introduced. 

This can be explained by the fact that the sediment feed was heavy and a con- 


94. “Elementary Mechanics of Fluids” by Hunter Rouse, John Wiley & Sons, 


(46) 


TABLE I 
Summary of Experimental Results of Properties of Sand Waves 
Run No. 2-1 2=2 2-3 2-4 2-5 
discharge cfs 0.4 0.4 0.4 0.4 
depth ft 0.174 0.190 0.215 0.236 0.267 
velocity fps 2-3 2.1 1.86 1.70 1.50 
slope 0.0115 0.0086 0.0054 0.0032 0.002 
Ry ft 0.035 0.051 0.06 0.077 0.092 
Rp ft 0.162 0.171 0.189 0.189 0,218 
Fy 0.972 0.848 0.707 0.615 0.512 
amplitude ft 0.089 0.078 0.059 0.037 0.053 
wave length ft 0.92 0.835 0.834 0.715 1.58 
wave vel. fps 0.0455 0.0236 0.014 0.0093 0.0051 
t° F 72 77 73 76 75 
Ty #/ft 0.116 0.092 0.064 0.038 0.027 
Vy fps 0.245 0.22 0.18 0.14 0.12 
Sedi. supply #/hr 600 300 150 75 38.3 
Sedi. Comput. #/hr 440 297 157 7265 = 
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ec, Velocity of sand waves in fps 


6 1 
V, Mean velocity of flow in fps 


Fig. 8 = Relationship between Velocity of Sand Wave and 
Velocity of Flow 
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Discussion of 
“MECHANICS OF STREAMS WITH MOVABLE 
BEDS OF FINE SAND” 


by Norman H. Brooks 
(Proc. Paper 668) 


T. BLENCH,! M. ASCE.—From the description of experiments the writer 
believes that a steady state of sediment movement had not been reached and 
that, even if it had, no quantitatively useful results could have been obtained 
with the suspended and bed loads measured as one. The test for equilibrium 
(“The two sets of profiles were practically identical, indicating that the equi- 
librium achieved was stable, well defined, and reproducible.”) seems to be no 
more than a test of the fact that the same conditions applied for the same time 
produce the same degree of disequilibrium. It is not surprising therefore that 
“the data obtained are not sufficiently extensive to permit establishing quanti- 
tative formulas.” Why they should be considered sufficient to destroy two 
well-established theories, of which regime theory is one, is not clear. 

The writer draws attention to the following misunderstandings of regime | 
theory. First that it is his own, whereas he has always been at pains to attrib- 
ute it to Gerald Lacey. Second that it avoids consideration of shear, whereas, 
in the way he presents it, gdS is brought explicitly into the regime slope for- 
nfula. Third that basic data were not made available, whereas the quoted ref- 
erence (15) gives its own references 2, 3, 6, 7, 9, 10 abounding with data, 
Fourth that it is based on an erroneous assumption concerning suspended load 
(vide SUMMARY OF CONCLUSIONS) whereas it makes no assumption about 
suspended load and offers no means for dealing with it; it is essentially a bed- “ t 
load theory. Actually the theory consists of explaining observed facts of canal Ke 
self-adjustment that were found to be represented by certain formulas; the 
formulas do not depend on the theory. 


1, Prof. of Civ. Eng., University of Alberta, Edmonton, Alberta, Canada, and | ' 
Con. Engr. 
748-47 
| 

| 


| 
} 
| 
| 
| 
{ 
4 


4 
= 
iz 
a> 


PROCEEDINGS PAPERS 


The technical papers published in the past year are presented below. Technical-division 
sponsorship is indicated by an abbreviation at the end of each Paper Number, the symbols 
referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics 
(EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), Sanitary 
Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping 
(SU), and Waterways (WW) divisions. For titles and order coupons, refer to the appropriate 
issue of “Civil Engineering” or write for a cumulative price list. 


VOLUME 80 (1954) 
JULY: 457(AT), 458(AT), 459(AT)°, 460(IR), 461(IR), 462(IR), 463(1R)°, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 
475(SM), 476(SM), 477(SM), 478(SM)°, 479(HY)°, 480(ST)©, 481(SA)©, 482(HY), 483(HY). 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)©, 488(ST)°, 489(HY), 490(HY), 491(HY), 
492(SA), 493(SA), 494(SA), 495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 500(HW), 501(HW)°, 
502(Ww), 503(WW), 504(wWw)°, 505(CO), 506(CO)*, 5O7(CP), 508(CP), 509(CP), 510(CP), 
S11(CP). 


OCTOBER: 512(SM), 513(SM), 514(SM), 515(SM), 516(SM), 517(PO), 518(SM)°, 519(IR), 520(1R), 
521(IR), 522(R)°, 523(AT)°, 524(SU), 525(SU)®, 526(EM), 527(EM), 528(EM), 529(EM), 
530(EM)°, 531(EM), 532(EM)-, 533(PO). 


NOVEMBER: 534(HY), 535(HY), 536(HY), 537(HY), 538(HY)°, 539(ST), 540(ST), 541(ST), 542(ST), 
543(ST), 544(ST), 545(SA), 546(SA), 547(SA), 548(SM), 549(SM), 550(SM), 551(SM), 552(SA), 
553(SM)°, 554(SA), 555(SA), 556(SA), 557(SA). 


DECEMBER: 558(ST), 559(ST), 560(ST), 561(ST), 562(ST), 563(ST)°, 564(HY), 565(HY), 566(HY), 
567(HY), 568(HY)°, 569(SM), 570(SM), 571(SM), 572(SM)°, 573(SM)°, 574(SU), 575(SU), 576(SU), 
577(SU), 578(HY), 579(ST), 580(SU), 581(SU), 582(Index). 
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JANUARY: 583(ST), 584(ST), 585(ST), 586(ST), 587(ST), 588(ST), 589(ST)©, 590(SA), 591(SA), 
592(SA), 593(SA), 594(SA), 595(SA)°, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(WwW), 
616(WW), 617(IR), 618(IR), 619(IR), 620(IR), 621(IR)©, 622(IR), 623(IR), 624(HY)©, 625(HY), 
626(HY), 627(HY), 628(HY), 629(HY), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(PO), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)°, 642(SA), 
643(SA), 644(SA), 645(SA), 546(SA), 647(SA)°, 648(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)°, 655(SA), 656(SM)°, 657(SM)°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)©, 662(ST), 663(ST), 664(ST)©, 665(HY)©, 666(HY), 667(HY), 
668(HY), 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 
€78(HY). 


MAY: 679(ST), 680(ST), 681(ST), 682(ST)©, 683(ST), 684(ST), 685(SA), 686(SA), 687(SA), 688(SA), 
689(SA)°, 690(EM), 691(EM), 692(EM), 693(EM), 694(EM), 695(EM), 696(PO), 697(PO), 698(SA), 
699(PO)©, 700(PO), 701(ST)°. 


JUNE: 702(HW), 703(HW), 704(HW)°, 705(IR), 706(IR), 707(IR), 708(IR), 709(HY)°, 710(CP), 
T11(CP), 712(CP), 713(CP)°, 714(HY), 715(HY), 71€(HY), 717(HY), 718(SM)°, 719(HY)¢, 
720(AT), 721(AT), 722(SU), 723(WW), 724(WW), 725(WW), 726(WW)°, 727(WW), 728(IR), 
729(IR), 730(SU)°, 731(SU). 


JULY: 732(ST), 733(ST), 734(ST), 735(ST), 736(ST), 737(PO), 738(PO), 739(PO), 740(PO), 
741(PO), 742(PO), 743(HY), 744(HY), 745(HY), 746(HY), 747(HY), 748(HY)°, 749(SA), 750(SA), 
751(SA), 752(SA)°, 753(SM), 754(SM), 755(SM), 756(SM), 757(SM), 758(CO)°, 759(SM)°, 
760(WW)°, 


c. Discussion of several papers, grouped by Divisions. 
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